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^ 1 Introduction 

Under rather general conditions on the kernel K(x, y), a symmetric homogeneous function in x and 
y, the Cauchy problem for the Smoluchowski coagulation equation: 

_ I = «[/] (i.i) 

X! 

O , i rx roc 

Q[f] = 2/ K(x-y,y)f(x-y)f(y)dy- K(x,y) f(x) f(y) dy (1.2) 

J J 

£ f(0,x) = f in (x)eL\(R), (1.3) 



> 
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o 

> 

x 



has a global solution f(t,x) for all initial data fi n (x) such that Jq 00 x fi n (x)dx < 00. Moreover, this 
solution satisfies the same estimate for all t > 0. 



Equation (1.1), (1.2) describes the aggregation process of particles of mass x and y with prob- 
ability K(x,y), assuming that the distribution of particles are uncorrelated at all times. In this 
— context the quantity: 

I xf(t,x)dx (1.4) 
Jo 

represents the total mass of particles in the system. 



On the other hand, it is known that, when the kernel is of the form K(x,y) = x a yP + x^y^ 
with a > 0, > and a + (3 = A > 1, the solutions to the Cauchy problem for the Smoluchowski 
equation undergo the so called gelation phenomenon. This means that there exists a positive time 
T g < 00 such that, for all t < T g , 

OO POO 

xf(t,x)dx= / x fi n (x) dx (1-5) 
Jo 
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and for all t > T g , 

rco poo 

/ x f(t, x) dx < / xfi n (x)dx. (1-6) 
Jo Jo 

It is also known that when A < 1 gelation does not occur and mass is conserved for all time (cf. 
[6j [H]). Two interesting open questions are related with this phenomenon. One is to describe the 
solution / as t — > T g . We consider here the second one, which is to understand the behaviour of the 
solution / after the gelling time T g . 

Althought no general result is known, several partial results indicate that before the gelling 
time, at least for a large familiy of initial data, the solutions to the coagulation equation decay 
exponentially fast as x —* +00. The Smoluchowski equation (1.1) has a discrete counterpart: for 
which an explicit exact gelling solution was constructed in |13j for K(k,j) = kj. Such a solution 
decays exponentially fast before the gelling time, and as a power law after that time. The exponential 
decay, before the gelling time, was later shown in [7] for the continuous equation (1.1), A = 2 and 
several initial data. Moreover, it has also been formally shown in [3] [7J that, for several initial data 



and A E (1)2], the solution of (1.1) decays, after gelling, like x ~( + >> as x — > +00 (see [12] for 
more detailed references). On the other hand, it was proved in [8] that £-( 3 + A )/ 2 i s the only possible 
power law decay for the solutions of (1.1) after gelation. Our main purpose is to prove that for the 
coagulation kernel 



K(x,y) = {xy) x '\ AG (1,2) 
and any initial data fi n , regular near the origin and such that: 

fin(x) ~ x~ (3+A)/2 as X -> +00, 



(1.7) 



(I- 



the problem (1.1)- (1.3) has a solution / satisfying 

f(t,x) ~ a(t) x~ (3+A)/2 as x +00. 



Moreover, this solution satisfies 



d 



dt 



x f(x, t)dx = -27ra 2 (i) for all t > 



(1.9) 



(1.10) 



which was formally shown in [5] for the discrete equation. By (1.10) the total mass of the solution 
/ is decreasing. This loss of mass is a characteristic feature of the solutions of (1.1), (1.2) after the 
gelation time. The choice of exponents A < 2 is natural, because A < 2 excludes instantaneous 
gelation or non existence of solutions ( [31 |5j [15]) and A = 2 is one of the "explicit" cases which has 
been treated using the Laplace transform (cf . [7J ) . 



In order to prove the existence of classical solutions of (1.1 )-( 1.3 ) after gelation we will use the 
same approach as in [2] [TO]- The starting point of this approach is to linearize around an initial 



data f in satisfying / if 



A )/ 2 f or x i ar g e anc j to derive detailed estimates on the solutions of 



the resulting linear equation. 



dg_ 
dt 



^{9 if in)- 



(1.11) 



To this end we will need some rather delicate estimates on the asymptotics of the solutions as x 
tends to infinity. Moreover, even to prove solvability of the linearized problem ( |1.11[ ) is nontrivial. 
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We will obtain it treating this problem as a perturbation of the problem obtained replacing f\ n by 
its asymptotics as x tends to infinity: 

% = m- <i.«> 



In order to carry on this program we need to derive detailed estimates about the solutions of (1.12). 
This will be the main goal of this paper. 

The linearised equation around the weak solution ;c~( 3 + A )/ 2 may be introduced more directly 
as follows. Consider a solution f(t, x) of the coagulation equation with an initial data fi n sat- 



isfying (1.8). If one is interested in the behaviour of f(t,x) for x large it is natural to scale 



the variables as follows: x = Rx,y = Ry,t = RT^'^/H and f(t,x) = i?-( 3+A )/ 2 F R (t,x). In 



these new variables, the equation (1.1) reads (-Fr)j = Q[-F!r] and the initial data Fj„ satisfies now: 
F R (0,x) = #( 3+A )/ 2 fi n (Rx) ~ (x)^+ x y 2 as R -> +oo. The limit of the function F R as R -> +oo, 



if it exists, would then solve the same equation ( 1.1 ) with initial data x ( 3 + A )/ 2 . Therefore the linear 



problem (1.12) appears naturally as the linearisation of the coagulation equation (1.1) in the region 
x » 1. Notice, however that in the region where x is small the function / is bounded and the 
approximation by means of the power law x ~( 3 + A )/ 2 cannot be valid. The analysis of that region 
would lead naturally to the study of a boundary layer whose description requires the analysis of the 
operator C. This will be made in a forthcoming work. 



On the other hand, the linearised equation ( 1.12 ) has some interest by itself. It is indeed a simple 



model to describe a set of particles at equilibrium, whose density distribution is given by a; - ( 3+A )/ 2 , 
and where a small set of particles is introduced, whose distribution <p(x) is considered as a small 
perturbation. The particles so introduced start to collide both between themselves and with the 
particles in the background. The equilibrium density distribution 

x -(3+A)/2 ig 

then perturbed. The 

distribution density function of the resulting set of particles may then be seen at any time t as the 



equilibrium distribution x~( 3+A )/ 2 and a remaining perturbation ip(t,x). The linear equation (1.12) 
only takes into account the collisions of the "particles in the perturbation" with the background and 
describes how the distribution of these particles evolves in time. It neglects the collisions between 
particles in the perturbation. This could be a reasonable approximation as long as the perturbation 
<p(t, x) remains small. Notice that the number of clusters in the background as well as the number of 
particles (the total mass) are infinite (since nor a;-( 3 + A )/ 2 nQ r X 1 ~( 3 + A )/ 2 are integrable in (0, +oo)), 
but the number of clusters and particles in the initial perturbation are finite. Our results show the 
following: 

• There is instantaneously an infinite number of "perturbed clusters", although their mass is 
finite. 

• As t —* +oo, the number of perturbed particles (the mass in the perturbation) tends to zero, 
but the number of perturbed clusters remains infinite. 

• The total flux of particles is perturbed at t finite but tends to the flux corresponding to the 
original equilibrium distribution as t — * +oo. 

Our results are obtained using classical Fourier analysis and the Wiener Hopf method, in a similar 
way as we did for the linearized Uehling Uhlenbeck operator in [9J although with an important 
difference. This is the regularising effect of the operator L, absent in the operator studied in [9], and 
coming from the fact that L is similar to the half derivative operator. The fundamental solution of 



(1.12) has then very different properties than that obtained in |9j. 
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In Section [2] we state our main results and transform the integro differential equation ( 1.12[ ) to 
a Carleman equation in the complex plane. In Section [3] we state the fundamental properties of 
the auxiliary function $ appearing in the Carleman equation. This equation is solved in Sections [4] 
and [5] using the classical Cauchy integral, which gives an explicit solution. In sections [6] and [7j the 
precise asymptotics of the solutions are obtained. The Section [8] is devoted to a brief mention of 
the initial value problem. Some properties of the fluxes of particles described by the solutions are 
considered in Section [9} We have finally added Appendices I, II and III where are collected some 
necessary technical results. 



2 The Linearized Equation 



We start this Section writing the precise expression of linearized equation (1.12). 



Proposition 2.1 The linearised equation of (1.1)-(1.2) with K(x,y) = (xy) A//2 around the solution 
G(x) = x-^ +x )l 2 is 



dg_ 
dt 

L{g) 



x/2 



(x - y) x/2 G{x -y)- x x / 2 G{x)) y x/2 g{y)dy 



(2.1) 
(2.2) 



+ 



x/2 



((x - y) x ' 2 g{x -y)- x x / 2 g{x)) y^' 2 dy 



_ x -3/2 / y^ g{y ) d y- 2 V2x^- 1 ^ 2 g(x). 

Jx/2 



□ 



Proof. By the symmetry of the kernel K: 

1 /-sb/2 

2 J K(y,x -y)f(y)f(x - y)dy 



K(y, x - y)f(y)f(x - y)dy 



(2.3) 



x/2 



Then we may then write the equation ( 1.1 )-( 1.2 ) as follows: 

i-x/2 



df 

dt 



(x - y) x/2 f(x -y)- x x / 2 f(x) y x/2 f(y)dy 



x/2 



K(x,y)f(x)f(y)dy. 



(2.4) 



If we linearize around the solution G(x) = x ( 3 + A )/ 2 ; define / = G + g and neglect quadratic terms 
on g we obtain (2.1), (2.2). 



Remark 2.2 The second term in the right hand side of (2.2) can be seen as some kind of half 
derivative operator applied to function x x / 2 g(x). This will appear again in the Fourier analysis that 
will be done later on the linearised equation. 



Remark 2.3 In order for the first integral in the right hand side of (2.2) to be defined we need 
y 1+x ^ 2 g(y) to be integrable at the origin. For the second integral we need some kind of regularity of 
g(x) with respect to x. For example y x ^ 2 g(y) ^-Holder continuous with 7 > 1/2. Finally, for the 
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last one we need y x ^ 2 g(y) to be integrable as y — > oo. Assuming power like behaviours we then need 
bounds on g of the form: 



g{y) < c y 



-A/2-r 



as y — > +00, 



g{y)<Cy- x l 2 -P as y^O, 

for some r > 1 and p < 2. 

Definition 2.4 We otZZ denote as V(r,p) the set of functions 

V(r, p) = { 5 G L fo c ( M+ ); SU P 5(a;)a; A/2+p < 00, sup <?(x)x A/2+r < 00 



(2.5) 
(2.6) 



(2.7) 



We state now the main results of this paper. The first one is an existence and uniqueness result of 
fundamental solutions for the equation (2.1), (2.2). 



Theorem 2.5 For all xq > 0, there exists a unique solution g(t,-,xo) of (2.1), (2.2) with initial 
data: 

g(0, x, x ) = S(x - x ) (2.8) 
such that g(t, •, xq) G V(3/2, (3 — A)/2) for all t > 0. Moreover, g(t,-,xo) has the self similar form 



1 

x ' 



g(t, x, x ) = — g (tx . 



2 -,11. 

x 



(2.9) 



The function g(t, ■, 1) satisfies the following behaviours in the different regions of the k,t space for 
some explicitly known constants a{, i = 1, ■ • ■ , 5 and arbitrarily small positive constants e and S. 



We have the following representation formula for t > 1: 

2 

g(t,x,l) = t^tpi(a) + ip 2 {t,a) 
where a is the self similar variable: 

2 

a = fA=i x , 

and the functions (pi and (f2 satisfy the following estimates: 



(2.10) 



(2.11) 



ai a 2 + O [a 



for < a < 1 



(2.12) 



{ a 2 a-^ +0 (a^ l+x+£ )) for a > 1 
where a\ and a 2 are two explicit constants and for some positive constant e arbitrarily small, 

' bi(t) a~l + (t^- 5 ° a~l +£ ) for < a < 1 



ip 2 (t,a) = < 



b 2 {t)a~^ +0 (t^- 5 ° 



(2.13) 



3+A 

a ' 2 



for a > 1 
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2 r 

where bi andb2 are two continuous functions such that \b\(t)\ + \b2(t)\ <Ct x ~ l , e > and 5$ > 0. 
On the other hand, for all < t < 1 fixed: 



g{t,x,l) = < 



a 3 £2r§ +6 3 (t)x~i +C (ix"5 +<51 ) for < x < \ 



(2.14) 



3+A 3+A / 3+A r . 

a 4 ix 2~+6 4 (t)x 2-+C 2 <M /orx>§, 



where 63 and 64 are continuous functions such that | £>3 (t)| + |o 4 (i)| < Ct 1+S , with 5 > 0, o~i > and 
03, a 4 explicit numerical constants. 
Finally, if t — > 0: 



' t- 2 * (^i) + (r 2 ) /or x = 1 + 0(t 2 



#(t,x,l) = < 
where the function ^ is given by: 



O 



for t 2 < \x - 1| < \ 



(2.15) 



§e * 3 ' 2 , / or a /z x > 0, 



(2.16) 



for all x < 0. 

Remark 2.6 For t > 1 /arc/e t/ie behaviour if g in terms of the variable x is : 

a 1 1^ cr~l + oi(i)cr~i + O (t*=i cr -1 ^^ + cr~§) /or < o- < 1 



#(i,x, 1) = < 



a 2 t^i + 6 2 (t)o--^ + O (t& o-( 1+A+£ ) + t^~ S ° o--^) /or a > 1 

(2.17) 



a x t ^x~%+0(t ^-ix~5 +e + t A-i d x f) for0<x<2t 



2 

A-l 



. 3+A 



A+1+2e 3+A 



^4-<5 -I 



(2.18) 



a 2 t A ~ 1 x 2 + O It a-i x 2 e + t a-i ° x 2 for x>2t 



Our second result is about the regularity of the fundamental solutions of (2.1), (2.2). 



Remark 2.7 Using the rescaling properties of the problem (2.1), (2.2), (2.8) has the self similar 



form (2.9). It is then enough to restrict our analysis to the case xq = 1. 



Our strategy to solve the problem (2.1), (2.2), (2.8) is to use Fourier analysis. The resulting 



problem is explicitly solvable by means of the Wiener Hopf method [2] . Using the representation 



formula for the solution, we then prove Theorem 2.5 by deriving suitable a priori estimates. Similar 
arguments have been used in [9] . 
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2.1 Fourier variables. 



We make now a change of variables in order to have functions defined in all of the real line M. To 
this end we define x = e , X G R, as well as the Fourier transform 



,6 = 4=/ e-^G(t, G(t, X) = e x ) 

V27T Jm 



(2.19) 



Then, the problem (2.1), (2.2), (2.8) reads in terms of the new variables: 



9G, ^ ~f , A-l\ / A-l. 



G(0,0 = 
where the function <]? is given by: 



i 



2tt 



2y^F r(z£ + 1 + 1; 



(2.20) 
(2.21) 

(2.22) 



as it is shown in Section 10 The fact that the function g(t, •, 1) G V(3/2, (3 + A)/2) implies that the 
function G(f, •) is analytic in the strip: S = {£ G C; Im£ G (3/2, (3 + A)/2)}. 

3 The auxiliary function. 



The properties of the function <£, defined by (2.22) determine the properties of the equation (2.20). 
The set of its zeros is: 



1 + A , 

i I n H — I , n = 0, 1, 



The singularities of the function are given by the family of poles: 

£ p (n) = i (l + ^ + n\ , n = 0,l,--- 

The asymptotic behaviour of the function $ as i2e(£) — > =too is the following: 
Proposition 3.1 For all M > fixed: 



(3.1) 



(3.2) 



*(0 



1 



as -Re(£) — > oo, uniformly on Im{£) G (— M, M) and where the function Q is defined as: 

Q = Q(0 = sgn(Re[^)). 



(3.3) 



7 



6 + A 



5 + A 



4 + A 



2 X V 


T , 3 + 3 A 
Im S= 4 


2 + A . . 
2 1 * 


> 


1 + A , M 
2 < 


J Im f= 4 


-7.5 -5 -2.5 

• Poles 

Q Zeros 


2.5 5 7.5 



Figure 1: Some relevant zeros and poles of the function <£. 



4 Solving (2.20)-(2.21) 



Our goal is to solve the problem ( 2.20 )-( 2.21 ) which we recall here: 

9G , - / , A - 1 .\ A A - 1 . 



G(o+,e) 



i 



2?r 



(4.1) 
(4.2) 



where the function G must be analytic in the strip S. This analyticity condition is a consequence 
of the decay required for G (t, X) in order for the integral equation to make sense. Since we are 
interested in deriving a solution G (t, X) in the sense of distributions, we want to obtain boundedness 
of G as |i?e(£)| — ► oo. In this case we will obtain, for the particular solution G constructed here, 
exponential decay, something that means that G(t, •) 6 C°° for t > 0. 



The following transformation allows to reduce (4.1 ), (4.2 ) to an equation for a function depending 
only in one variable £ : 

V2 r 

Ilm(y)=f3 V (y) ' 



G(t,0 



V (£) 2»u- W ) 



A- 1 



/tt* (A — 1) 
for some (3q G (3/2,2). We have: 

Lemma 4.1 Suppose that V {rj) is analytic in the strip S and that V satisfies 

e~V^^ [ y/\y\ + l\dy\ < oo 



dy 



(4.3) 





1 


/ 

Um(y)=Pa 


V(y) 



(4.4) 
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for any /3q G (1)2), as well as: 



V(i/) = -V U + 



A-l 



A-l 



(4.5) 



for Im{rj) G (§,2) . Define G (t,£ ) by m eans of (4-3) for im (£) > (3$. Then G can be extended 
analytically to S and it solves (4-D, (4-%) f or Irnirj) G (§,2) . 

Proof. It just follows by direct computation. Indeed, notice that Stirling's formula, that is 
uniformly valid for V {z) , arg (z) G (— ir + Eq, tt — Eq) with £q > (cf. [1]) implies: 



2j (e - y) 
A-l 



e ix=iji"i 



for |£| < R, Im{y) = j3$. Therefore, th e int egral on the right-hand side of (4.3) converges for any 
^ G Sn {£ : Im ({) G (/3 , ^r) } due to (O and the function G {t, £) satisfies: 



G(t,0 



<C R 



Im{y)=Po 



V(y) 



c -rr^\y\ \ d y\ 



< oo , \£\<R 



Taking (3q arbitrarily close to | we obtain analyticity of G in S. 



Moreover, the derivative with respect to t of G in (4.3) can be computed by means of: 

V2 



dG 



(t,0 



V(0 . 

TTi(X - 1) J Im(y)=f3o V(y) 



A-l 



-1 



2j (g - y) 
A-l 



+ 1 My 



(4.6) 



where we have used zT {z) = T (z+ 1). On the other hand, using (4.3) we obtain: 



(A-l) 



7U (A — 1) Jj m ^ = p 



Mi=v) i,, / 2i($-y) . , 
* r ( ^ _ ^ + 1 ) 



and using (4.5) and (4.6), (4.1) follows. 



It only remains to check (4.2). To this end we use contour deformation and residue Theorem to 

G r (t,0 



transform (4.3) into: 

G(t,0 



1 



2tt 



G r (t,0 



v/2 



V^i (A - 1) Ji m (y)=p 1 V (y) 



t r 



A-l 



dy 



where (3\ > Im(^) . Using (4.4) it follows that: 

2Q3 1 -7m(£)) 



G r (t,£) 



<c R t 



£GSn{|£| <R} 



and therefore G r converges to zero uniformly in bounded sets of £, whence (4.2) follows. 



□ 



Remark 4.2 A heuristic explanation for the formula (4-3) can be given using Laplace transform. 
Suppose that we define the Laplace transform of G(t,£) in t as: 



G(z,0 



G(t,£) e~ zt dt 
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Then, (4-D, (4-%) becomes: 



A-l 



zG(z,0 = G t,e+— -i *U + — -* + 



A-l 



1 



/2vr 



(4.7) 



TTie solution of this equation can be formally reduced to (4-5) by means of the transformation: 

2i 



G(z,g) = exp 



A-l 



(4. 



The reason for using log (—z) instead of log (z) is that can 6e expected to be analytic for 

Re (z) > a for some o£l and in this form the function log (—z) can be expected to be analytic in 
this region assuming that arg(z) G (— tt, 7r) . The transformation (4-8) brings (4-D to: 



H(z,£)-Hlz,S+^—±i 



1 e j*r ]<*(-*)* 
~2^ zV (£) 



(4.9) 



Equation (4-9) can be solved using conformal mapping and Cauchy's formula. Using the change 
of variables: 

H(z,£) = h{z,() , C = e^~ M (4.10) 
and where, for the sake of simplicity we will write, with some slight abuse of notation 

v(0 = v(0 

we obtain: 

h{z,C + i0)-h(z,(-i0) = = \— , ( £ R + (4.11) 



f%H zV (C) 



with h analytic in C \ M. + and: 



a (z) 



2-iri 



log (~z) 



mates for small and 



The solution of (4-H) can be obtained, assuming that y^y satisfies suitable boundedness esti- 



Large Q, and using Cauchy's formula: 

. s 1 1 e^ oa{z) f°° s a ^ ds 
h(z,Q 



and, using (4-10): 



H(z,0 



2m ^2k : J V (s) (s - Q 

dy 



1 1 1 

2iri J2tt z 



Aira. (z) 



(4.12) 



It then follows from (4-8) that: 



G(z,0 



1 1 1 

2ni y/2~TT z 



v(0 



4:7YOl(z) / /-\ 



dy 



V(y) 



1 - e A 



and inverting the Laplace transform we finally obtain (4-3). 



10 



5 On the solutions of (4.5). 



Equation (4.5) admits infinitely many solutions. Indeed, given any solution V par t (£) we can obtain 
any other one by means of: 

V (0= Vpart (Z)P(0 



where: 



P(Q = PlZ+t—±i 



(5.1) 



Given such a non uniqueness a natural and essential question is then how to chose one of them. 
We may state several sufficient conditions that would ensure that G is the Fourier transform of 
a tempere d dis tribution. First we want the function G to be defined. This is guaranteed by the 
condition (4.4) above. However, this condition is not sufficient to prove that G(i, £) is globally 
bounded with respect to £. The difficulty comes from the fact that, if the behaviours of V(£) are too 



disparate as Re((,) tends to plus or minus infinity, the quotient 



v(0 
v(?+y) 



may be strongly increasing 



in some regions of the integral in (4.3). A sufficient condition to avoid this difficulty is to have: 

TT 



\v(0\ 



\B±\ < 



(A-l) 



(5.2) 



as Re(£) 

growth of the quotient 



±oo. The decay rate of the Gamma function in (4.3) may then control the possible 
v(0 



uniformly on £. 



v(5+r) 

Another requirement that we need for the function V comes from the requirement that G must 
be analytic in the strip S. This is ensured by imposing also that V is also analytic in that same 
strip. 

We will then construct a function V analytic in that strip, satisfying equation (4.5) for Im(£) 6 
(3/2,2), satisfying conditions (4.4) and (5.2). 

A particular solution of ( 4.5[) ca n be easily obtained using Cauchy's formula. To this end we take 
the logarithm of both sides of (4.5) to obtain: 



log(V (£)) = log( V(£ + ^-i 



+ log 



(5.3) 



or equivalently 



iog(v (0) + log (-$(£)) • 



(5.4) 



Let us take any (3\ such that $ (£) has no zeros nor poles along the line Im (£) = 0%. We define: 
^(C) = log(V(£)) , C = e ^-^ , Q(C) = log(-*(e)) 



Equation (5.4) then becomes 

iP(C + iO) =^(C-iO) + Q(C-*0) , CG 



(5.5) 



with ip analytic in C\M + . Taking into account that |Q(C)| — C (1 + |log (C) I ) we can obtain a 
particular solution of (5.5) as: 



<M0 



i 

2vri 



Q(s) 



l 



l 



ds 
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where the term l/(s + 1) has been added to the classical Cauchy integral in order to ensure the 
convergence of the integral. Then, returning to the variable £ we obtain: 



V, 



part,{3i 



(£) = exp 



(A — 1) i Ji m (r))=p : 



log (-$(7/)) 



1 + e 



dq 



(5.6) 



where Im (£) G (ft 



A-l 

2 • 



ft) 



Formula (|5.6|) provides a particular solution of (4.5 ). On the other hand, we can obtain an infinite 



family of solutions of (5.1 ) given by: 



Pit) 



4-7T , 

e A-i* 



Let us define a family of solutions of ( |4.5[ ) : 

V (0 = e^% arttl3l (£) 



(5.7) 



(5. 



Actually, using Fourier series, it can be seen that any solution of (5.1 ) can be written as a infinite 
linear combination of the functions Vg(£). 



The formula (5.6) does not define uniquely the function V pa rt,/3i unless we prescribe the value of 



ft and the argument of the function ln(— $(77)). The different possible choices of this argument just 
differ by a factor 2irl i and therefore the resulting functions V par t $ x would differ by a multiplicative 
factor (5.7). Proposition 3.1 implies that arg(— $(77)) — > ir/4 + 2ir£i as Re{rj) — > +00. In order to 
avoid exponential factors in some of the forthcoming formulas, we determine uniquely the function 
ln(— 3>(r?)) by choosing: 



lim arg( — $(77)) 



7T 

4 



(5.9) 



Notice that in the formula (5.6) there exists an infinite possibility of choices of the constant 
01. These functions may be extended analytically moving £ and simultaneously the contour of 
integration in such a way that the condition —(A — l)/2 < Im(£ — rj) < always holds. The only 
true obstruction to extend analytically the functions V par t ij g 1 arises from crossing with the contour 
deformation the zeros or poles of the function <]?. Suppose that tsing 1S a zero or a pole of $ and ft, 
ft are such that: 



sing 



- < ft < £smg < ft < £sing + - • 



Then: 



V, 



part, (3 



i(0 



where: 



47T ^ 

1 + e^-i^" 



1 if £sin 9 is a zero 
-1 if tising is a pole. 



(5.10) 



(5.11) 



Combining (5.10) with (5.6) we can then extend any function V par t,f3 to the whole complex plane as 
a meromorphic function. As it could be expected the different functions V par t,f3 can be related to 
each other by means of linear combinations of functions of the form given in (5.8). 
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In order to obtain the function V(£) with the properties requested above, it is sufficient to take 

A — 1 



V(0 = (0, M£)e A 



: ,ft 



with 



ft € 



2+A 3+A 



(5.12) 



(5.13) 



Moving the contour of integration if needed, inside the strip Im{rf) G (2 + A)/2, (3 + A)/2 we obtain 
that the function V(£) has no zeros nor poles in the whole strip S. The function V is then defined 
once and for all along the remainder of the paper. It only remains to check that this function satisfies 



the two conditions (4.4) and (5.2). It follows from Proposition 11.2 in Appendix 11 that: 

Cse -U^m < | V(f) | < Cse -k(xh-m\ for Jm(e) e f|,i± v 



(5.14) 



for 5 > arbitrarily small and C$ > a constant depending on 5. This behaviour implies both (4.4) 



and (5.2). 



Summarizing, we have shown: 



Proposition 5.1 The function V(£) defined by means of (5.6), (5.12) with ft as in (5.13) can be 
extended analytically to the strip S and meromorphically to the whole complex plane. It satisfies the 



equation (4-5) as well as the estimates (5.14). Morever, V(£) ^ in all the strip S. 



Proposition 5.2 The function G(t,£) defined by means of (4-3) with V defined in Proposition 



5.1 



solves (4-1) (4-%)- Moreover we have the following representation formula: 

In (— <J>(r?)) x 



v(0 
v(iO 



exp 



(A - l)i J Imr)= p 1 



1 



1 



dr] 



(5.15) 



for £ and y such that Pi - (A - l)/2 < Im(^) < 0i and 0i - (A — l)/2 < Jm(y) < ft. 

6 Decay estimates for the function G(t,£). 

Through all the following of this paper we shall use the following function A{z) that we define by 
means of: 



T(z) = V2^e- z z z ~ 1/2 A(z). (6.1) 

Notice that by the Stirling's formula: 

A(z) -» 1 (6.2) 
uniformly as \z\ — > oo and arg(z) G (— tt + Eq, tt — Eq) for any £q > small. 
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6.1 The case < t < 1. 



Lemma 6.1 The function G defined in (4-3) satisfies the following. 

\G{t,Q\ < Ke' a ^ 1 



i + iei 



1/2 



d_ 



G(t,0 



+ u + iei 



3/2 



< K te- a Vl?l* 



(6.3) 
(6.4) 



for all t G (0,1) and for all £ G {£ G C; Im(£) G (3/2, (3+A)/2)}U{£ G C; Jmf G [-L, L] , |f| > 2 L} 
for any L > sufficiently large and where K\ and a are positive numerical numbers depending on L. 

The proof of this Lemma requires the following result: 
Lemma 6.2 For any fixed constant B > consider a function 

i(Z,Y)e*& Y >V dY 



W(t,0= [ 



ml 



' Irn(Y)=—yi 

where the function m(£, Y) is analytic with respect to Y in the domain: 



(6.5) 



Re I Y + -sign (-Re(£)) £ 



< B\Im(Y)\, sign(Re(Y)) = sign(Re(£)) 



and in the strip: 



ImiY) G 



-7i ,7i + 



A-l 



Then, for any L > sufficiently large, there exists £o > sufficiently large, and depending on L, 
such that: 



\m{Z,Y)\ < C 



then, 



\W{t,0\<Ce 



-a\H\ 1/2 t 



(6.6) 



(6.7) 



for all t G [0, 1] and all £ such that \Re(£)\ > £ and {£ G C; Im(£) G (3/2, (3 + A)/2)} U {£ G 
C; ImCe[-L,L], \£\ >2L}. 



then, 



\m(£,Y)\ < C(1 + \Y\) and m(£,0) = 



|W(t,0| < Ct\i\ l ' 2 e- a ^ 1/H 



(6.8) 



(6.9) 



/or a// 1 G [0, 1] and all £ suc/i t/iat |J2e(£)| > £ and {£ G C; Im(£) G (3/2, (3 + A)/2)} U {£ G 
C; Im£ G [-L, L] , |£| > 2L}. 
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ImZ 



ImZ = 7x t 



ReZ 



Z c 



Proof of Lemma 



6.2 



function ^ becomes: 

' 2iZ 



Figure 2: The curve C% 
It is convenient to introduce the new variable Z as: Y = y/\£\Z. Then the 
2 



2iZ 

131 ~MH y 



(A-l)i 
2iZ / 2iZ 



In (-<%)) 6(77 -^y^Z)^ 



a-i^ + a^i-Ia^i"^ 



1 



ln|^| 1/2 



(6.10) 



The possibility of extending the function Z, t) as a function of Z to some suitable cones has 
been studied in Lemma 12.1 Moreover in Lemma |12.3 the existence of a critical point Z c of the 
function Z, t) with he asymptotics (12.12) has been proved. 

Using the analyticity properties of the functions <£(£, Z, t) and A{z) as well as the change of 
var iables Y = y/\£\Z we can obtain a new representation formula of G 



W(t,0 = VW\ [ e^'^mfo VW\Z)dZ 
Jc t 



(6.11) 



with the new integration contour Ct given in Figure 2. Notice that we have moved the portion of 
the contour ImZ = (/?o — -f£)|£| _1//2 where |i2e(Z)| is large to the line Im(Z) = 71 since in that 
deformation we do not cross any singularity of the integrand. 

We now consider separately the two different cases |£| 2 t — ► +00 and |£| 2 t bounded. 



The estimate of W(t, £) as |£| t 2 — > +00. We write the function W(t,£) as follows: 

w(t,o =ii+i 2 
ii = 



12 



III / e*& z tim(£,VW\Z)dZ 

C t n{Z;\Z-Z c \<8 t} 

\i\ [ e^ z ^m^y^\Z)dZ 

C t n{Z;\Z-Z c \>8 t} 



(6.12) 
(6.13) 

(6.14) 
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for some positive constant 5q to be fixed. 



We estimate first the integral X\ using Lemmas 1 1 2 . 1| - [1231 and Taylor's expansion we obtain first: 

lei 



m Z, t) = m Z c , t) + {Z / c)2 S(^, Z c t) + O { ^i(Z Z, .) :! | (0. 1.5 ) 



dz 2 



for |Z — Z c | < egt, and then 



[l + 6(C,Z,t))\Z-Z, 



2 
C| ) 



(6.16) 



2 v^7Tt(l +iQ) 

where |<5(£, Z, t)| can be made arbitrarily small if |^|t 2 is large and Eq sufficiently small. Therefore: 



C t n{Z; \Z-Z c \<8 t} 



which gives, in the case (6.6) : 



Xi(t, 01 < Ce" a vl5l* as |^| t 2^ +00) o<t<l, 



and in the case (6.8) 



|2i(t,0l < Ct^e- a Vl?l* as |£| t 2 -» +oo, < t < 1. 
We must now estimate the integral I2 given by ( 6.14[ ). 



(6.17) 



(6.18) 



12 = VW\ [ e^ z ^m^,y/W\Z)dZ 
Jc t n{Z;\z~z c \>8 Q t} 



In order to estimate the integral I2 we use Lemma 12.8 To this end we split the integral as follows: 
I2 = Z21 + 2^2,2 



22,1 = VW\ I 

Jc 



c t n{z-, \z-z c \>s t}nj(M,t) 



m(Z,y/\£\Z)dZ 



2^2,2 

22,3 



lei / _ e^> z >Mz, VW\z)dz. 

C t \ 7 (M)n{Z; \Z\<e ly M} 



1^1 / _ e*«' Z '')m(e,v^Z)dZ 

C t \ 7 (M)n{Z; |Z|>eiv1ei} 



that, in the case (6.6): 



where j(M, t) is the portion of the curve Ct along which Im{Z) < 71 1. Using Lemma 12.8 it follows 

m(£, VW\Z)dZ 



< Ce^ z ^e- £Q VW\t^ f 

Jc t 



t n{Z; \Z-Z C \>S t}rr/(.M,t) 

< Ce^ z ^e- £ °VW\t^ t < Ce *(^ 01 t) e - £ o/2^t ; 



and in the case (6.8) : 



1 2A \ < ae^ z - t ) e - e °v / l«|t(^t)2<c( v 4^t) e ^ z -*) e - e o/ 2 v / ICI* 



(6.19) 



(6.20) 
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The estimate of X2 2 follows using Lemma 12.9 In the case (6.6) we obtain: 



X 2 , 2 | < C vUl f e~ a VW\\z\ dZ = Ce -aM^W\t_ 

J\Z\>Mt 



and in case (6.8) we have: 



X 2 , 2 | < C vUl f e-*/MW VW\ZdZ < C(^/W\t)e- a 

J\Z\>Mt 



(6.21) 



(6.22) 



since |£| 2 1 is large. 

The third integral 12,3 is estimated using Proposition 11.2 To this end we use the variable Y 
and the identity: 



r / 2»y 

1 1 A-l 



V(Z + Y) A (2iY 



where A is defined by formula (6.1 ). This reduces the estimate of I2 3 to the estimate of the integral: 



J 



/m(y)=71,|y|>ei|£| 



m(S,Y) 



v(0 



V(Y + 



_ 2iY 

t 



2iY 
A-l 



(6.23) 



Using that 71 > 0, < t < 1 and Stirling's formula, it follows that, in both cases (6.6) and (6.8): 

v(0 



J< / (1 + |F| 

'/m(K)=7i,|y|>ei|£| 



v(y + o 



e~^ |y| dY~. 



Proposition 11.2 gives the following bounds 
J C 



'/m(y)=71,|y|>ei|e| 
< C< e -2(X^T) e il€l e ^ICI. 



(6.24) 



with £q > if 5 is sufficiently small. If we are in case (6.8), we estimate: 



which is estimated by the right hand side of (6.9). 



This ends the proof of the estimate of I2 an d then of the Lemma 6.1 in the domain where |£| t 2 
large and < t < 1. 



The estimate of W(t,£) for |£| t 2 bounded. We suppose first that condition (6.6) holds. We 

(6.25) 



split the integral in (6.32) in two pieces: 



w(t,£) = j 1 + j 2 
Ji(t,0-- 



Irn 



Y^ lt ,\Y\<e \C\ Ate) + 



v 7 ^ 



A-l 



//< - I m(e,y) V(Q (2iY n 

2{ ' >/^(A-l)yj rn Y=- 7lt ,|y|> eo [ € [Ate)V(y + VA-1 
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The integral J2 is estimated with the same argument used to bound the integral J in (6.23). 

|J 2 | < Ce~ £0 ^ 

We rewrite J\ as: 



(6.26) 



where 



*(£,CM) 



lei* 



M0=— %.lfl<eo- 



2< 
A - 1 



1 



In 1 |£| 



1/2 



1 - In 
1 



2*C 
A — 1 

2 i C 



In . 
2 V A - 1 



+ In [2 ^ne lQ * 
+ /i(£,CM)- 



(6.27) 



Notice that Lemma 1 2 .2 1 implies 

IM£,CM)|<c f|d 2 t 2 + o 



1 

lei 

for t |C| < eoyl?!' e -^(£1 l£l — £0 and for some positive constants C = C(B) and £0 = Co(-B) 
sufficiently large. Notice that 

fte(*(*,CM)) < ^t^C- ^C^-lln (t\t\ 1 / 2 )+ Ce ^M\t\C\ 



< VW\t(c 



7T 



2(A-1) 



ICI 



1 



In t ICI 



for some positive constant C, for £ such that Im(Q 



71 



1/2 



and I C| < £o"^^> assuming that Eq is 



small. Using that ty |£| is bounded it follows that, in the case (6.6) 

\Ji(t,o\<c(VW\t) 1/2 [ 
v y j 1 

and computing the integral: 

\Ji(t,0\<c 



//m(C)=-^,|CI< £ 0^ 



Combining this estimate with (6.26), estimate ( |6.7| ) of the Lemma follows. 

In the case (6.8) we deform the contour integration in formula (6.5) using the analyticity properties 
of the function m(£, Y) as well as the fact that ra(£, 0) = to obtain: 



tt(A-I) 
A(-l) 

+ 



m ( £, — — 2 



V(0 . 



m 



/ m (y)= 7l + A_^ 



vr(A-l) 
A(-l) 



m I £ 



A- 1 



A- 1 



i + ch 



7m(y)=7i + ^ L 



2 y V 2 



(6.28) 
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using (4.5). By (6.8), m(£, — ^— i) is uniformly bounded. Using then Proposition 

|^i|<C(l + |£| 1/2 )i. 



3.1 



we obtain that: 
(6.29) 



We are then left with the integral K2 which is formally very similar to (6.5) although the integral 



contour is different. We split this integral in two terms J\ and J2 like in (6.25). The term J2 is 



bounded as (6.26) by a similar argument as before with one small difference which that the term 



m(£, Y) is now bounded as (1 + \Y\). The term J\ can be written as: 



1 in 



2iY 
A - 1 



dY. 



Using Lemma 



12.1 



^ |byCe(^) +5 °) |y| 



in the domain of integration. On 

2iY ,j. 



we may estimate | v(y+£) 

A — 1 Mr 1 I - rl 

the other hand, since ImY = 71 H — ^~ we may estimate t x - 1 by i A ~ 1 . Therefore, using the 
decay of the Gamma function: 



Jl(t,0\ < Ct 1+ ^ [ e^ +5o > l e-^d\Y\ 
J/ m (y)= 7l +^ 



(6.30) 



(6.31) 



the proof of the Lemma 6.2 



Combining (6.29) and (6.30) we deduce the estimate (6.9) when is bounded. This concludes 



□ 



Proof of Lemma 6.1 . Using the change of variables: y — £ = Y in (4.3) we obtain 

y/2 f 



G(t,£) 



where 71 is a positive constant sufficiently small. We rewrite the function G(t,£) as follows. 

V2 r 

^ri(\ - 1) Ji m (Y)=- 71 



G(t,0 



(6.32) 



(6.33) 



where 



(A — 1) i Jlmr)=fix 



ln(-$(?7))e(77-£,y)d?7 



2iY 
~A-1 

e(<7,r) 



ln(t) 



2^y / 2fy i\ / 2iy 

A^T + A^T-2 r A^l 



1 



1 



1 --r^jv 1 tt 2E t(-^+^) 

1 — e A -! 1 — e A -! k ' 



(6.34) 
(6.35) 



In order to obtain estimates (6.3) and (6.4) for bounded values of £ we use contour deformation. In 
particular, crossing the pole at Y = in integral (6.33) and using residue's Theorem we obtain: 

1 



Gi(t,0 



v(0 



v^*(a-i) y /my=71 v(y + o 



2»y 

t "f 



2^y 

A - 1 



dY. 



(6.36) 
(6.37) 
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with 71 > sufficiently small. Using Proposition 11.2 it follows that 

|Gi(t,0l<Cts3r 



uniformly for £ in bounded sets. This yields estimate (6.3) for £ in bounded sets. If we differentiate 

V2 



in (6.37) with respect to £ we obtain: 

—G(t,0 = -^ i{x 



1) / 



m V=7i 



m^,Y)t~^T 



2iY 
A - 1 



v(0 



5^ 



(6.38) 
(6.39) 



Using the analyticity properties of the functions Y) we deform the integration contour in 

the integral (6.38). The first singularity that is met is the pole of function T located at 

y = (A — l)ij2~. (This point is below the first zero of the function V(£ + y) which is located at 
y = (2 + (A/2) - f) % and 2 + (A/2) - £ > (A - l)/2.) We deduce 



^ G(t '°-^(A-1 
V2 



2V2^ _ A-l 



vri(A - 1) J ImY=11+ x=i 



m e (£,Y)r^T 



2iY 
A-l 



(6.40) 
(6.41) 



using the analyticity properties of the functions me(£,,Y) we deduce estimate (6.4) for bounded 
values of £. 

In order to complete the proof of Lemma 6.1 we apply the method of the stationary phase. To 

exp: 

V2 



this end we differentiate the expression in (6.33) with respect to £, obtain: 



G(t,0 



m (A - 1) J Im 



w 



(H,Y,t)A 



2iY 
A-l 



(6.42) 



for £ = 1,2 and use Lemma 6.2 with three choices for m(t, £): 

V2 



m(c,y) 
m(c,y) 

<20 



Tri(A-l) \X-1 



2iY 



V2 



iri(\ - 1) 



for £ = 1,2. (Notice tha t y i) and y i) are independent of t) 



5V, 



(lfl + i)^tt,y;t)ii 



2iy 



A-l 



(6.43) 
(6.44) 



The function m in (6.43) satisfies condition (6.6). On the other hand, using Lemma 12.10 



it 



follows that for the choices (6.44), i = 1,2: |m(£,y)| < C|y|. Moreover, by the definition of the 

dip 



function * (cf. (|6.34|» 



d_ 



In 



V(0 

v(e+y) 



and therefore m(£,0) = with the choice (6.44) and 1=1. This follows by a similar argument 
for the choice (6.44), 1 = 2. Applying then Lemma 6.2 the estimates (6.3) and (6.3) follow for |£| 
sufficiently large and this concludes the proof of Lemma 6.1 □ 
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ImZ 



ImZ = 71 



ReZ 



Z c 



Figure 3: The curve V 



6.2 The case t > 1. 



Lemma 6.3 The function G defined in (4-3) satisfies that for any Eq > arbitrarily small, there 
exists two positive constants k\ and a such that 

\G(t, 01 < Ki t~ ^-i +£o e- a VM ( 6 .45) 

for all £ such that Im(£) G (3/2, (3 + A)/2) and all t > 1. 



Proof. We deform the integration contour in the expression (6.32 ) to the line: Im(Y) = (A— 1) /2—eq 



with so arbitrarily small in order to avoid the zero of the function V(£ + Y) at £ + Y = 1. Then 
we perform the change of variables Y = and perform a new deformation of the contour of 

integration in the Z variable to the curve T> in Figure 3. 
to obtain: 



Then, 



G(t,Z) 



\G(t,z)\<Cy/\t;\t-—i 



a — 1)7© V ^ — ^ J 



iri (A 



(6.46) 



v 



MS,z,i) A [ 2 Vkl^ 



A-l 



We argue now in the same way as in the previous case with t = 1 splitting the integral in the same 
pieces to obtain the estimate of the Lemma. □ 



7 Estimates on the function G(t,X). 

We may now take the inverse Fourier transform of G(t, £) to obtain the function: 

' ' e iX ^G{t,t)dt 



V27T Ji m (£)=i3 
where (3q is in the interval (3/2, (3 + A)/2) 



(7.1) 
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Proposition 7.1 For all t > the function G(t,-) defined by (7.1) belongs to C°°(M) and, for 
I = 0, 1, 2, • • •. Moreover, for any fixed R > 0, it satisfies: 



d e G 



8X e 



(t,X) 



< 



R 



for < t < 1, |X| < i?. 



£2(1+*) 

and, /or any £o arbitrarily small and R > 0, i/iere exists a positive constant Ce te0i R suc/i i/iai 



ax 1 



(t,X) 



< ° e : £(hR . for t>l, \X\<R 



-£0 



(7.2) 



(7.3) 



Proof. This Proposition is just a consequence of (7.1 ) as well as from (6.3) for t G (0, 1) and (6.45) 
for t > 1. □ 
The estimates (7.2) and (7.3) do not give any detailed description of the function G. We now 



proceed to obtain detailed descriptions of the function G on the different regions in X and t. 
7.1 Behaviour as t -> 0, < \X\ < 1. 

We begin by showing that the function G behaves like a mollifier of the Dirac measure when t — > 
and for X small. The asymptotic profile of the solution G(t,X) can be heuristically guessed as 
follows. If we assume that the function g, solution of (|2.1[) with g(0, x) ~ 6(x — 1) is small, then 



equation (2.1 ) may be approximated, for t and x — 1 small, as 

" 1/2 g(x - y) - g(x) 



dg_ 
dt 



II 



3/2 



dy, g(x,0) = 5{x - l) 



(7.4) 



The equation in ( 7.4 ) describes the probability distribution for a system of particles whose size is 



initially equal to one and have a probability per unit of time Xo<y<i/2 ^3/2 °f coalescing with a 
particle of size in the interval (y, y + dy). If we write the function g as: 

/ X 1 

g{t,x) 

Applying Fourier transform we obtain: 

dg 



g{t,x)e ikx dk 



/2tt 



-i k 



dt 

m{k) 



(t,k) = m(k)g(t,k), 

1/2 e -iky _ -y 

n y 



g(k,0) 



2vr 



3/2 



dy. 



The multiplier m{k) can be computed explicitly but its exact formula is not needed to compute the 
asymptotics of the function g{t,x). The only relevant information that we really need is 



m(k) « — 2V vr k i, for \k\ — > +00. 



Then, 



g(t, k) 



-2v it kit 



/2vr 



for \k\ 



+00 



and inverting the Fourier transform we obtain 

g(t,x) 



t 2 V t 2 



as t — * 0, 
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with 



e x for x > 0) 



(7.5) 



for x < 0. 



Finally, since £ = e^-, £ — 1 = e x — 1 ~ X when X ~ we obtain 

as t -> 0. (7.6) 

The fact that the support of G(t, X) is contained in M + is a consequence of the interpretation of 
(7.4) in terms of coagulation of particles given above. 



7.1.1 The region where X = 0{t 2 ). 



We now prove that the fundamental solution G(t,X) behaves in the self similar form (7.6) in the 
region where X = 0(t 2 ) using the explicit representation formula given by (4.3), (5.15) and (7.1). 
We are interested in the limit lim f _ >0 + t 2 G(t, t 2 x) = ^(x) m compact sets of x- To this end we 
rewrite (7.1) as follows: 



G(t, X ) 



1 



2vrt 2 



Im{ V )=/3 t 2 



t 2 



drj. 



(7.7) 



Formula (7.7) suggests that in order to compute \P(x) we need to obtain lim t _ >0 + G if: f° r V such 
that Im(rjJ= Pot 2 . To this end we use the expression (6.32) that implies: 



G{t^/t 2 ) 
Proposition 7.2 



V2 



TTl I 



V(r?/t 2 



(A-l) 7 /m y=_ 7l V(Y + rj/t 2 ) 



2iY 

t 



2iY 
A-l 



dY. 



(7.8) 



lim ,(t 2 G{t,t 2 x)) =Hx) 



unformly for x in compact sets o/K where ip(x) * s as Jn ( 7-5). 



The proof of Proposition 7.2 requires several Lemmas. 

Lemma 7.3 For all Eq > and M > f/iere exists a function h £0) M{t) such that 



lim h £0tM {t) = 



and 



V(7?/t 2 ) 



V(F + r//t 2 ) 



_ 2i Y 

e * 



^ln(2 0Fv^) 



e ,M 



(t) 



(7: 



(7.10) 



for all Y such that \Im(Y)\ < 1/4, |i?e(Y)| < t ^ Qt ^ , for Im(r)/t 2 ) in compact subsets of (3/2, (3 + 
A)/2) and Eq < \rj\ < M. Moreover there is 5o > sufficiently small (depending on eq and M) such 
that: 

V( V /t 2 ) 



_ 2 i Y 

-t 



3tt \Y\ 

< Ce 4 a— i 



V(Y + r7/t 2 )" - v ' (7,11) 

for all Y such that \Im(Y)\ < 1/4, |-Re(Y*)| < fy, for Im(r]/t 2 ) in compact subsets of (3/2, (3 + A) / '2) 
and Eq < M < -W. 
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Remark 7.4 In the Lemmas 7.3 until 7.6 and in Proposition 7.2 we choose the branch of the 
function square root as follows: 



r z = I^VVi with e (- 



Proof. Lemma 



(5.15 ) we have: 



7.3 



V(r,/t 2 ) 



V(r]/t 2 + Y) 



is a consequence of Lemma 12.1 where £ = rj/t 2 , Z = Y/y/\£\. By formula 

ln(-$(ry))e(ry-r 7 /t 2 ,y)d?7 



exp 



exp 



and using (12.4): 



(A — 1) l Jlmrf^ 
F(jl/1?,tY/y/h\ 

F (r]/t 2 , t Y/ v^|) + In (2 ^7^? 



Therefore, in the region where | J?e(V)| < l/(t|lni|) we obtain (7.10). On the other hand, if 
\Re(Y)\ < 5 /t 2 we deduce (|7.1l| using that: 



Re 



2iY 



In I 2 V7ry i 77 



< 



7T 



2(A-1) 



A — 1 

as |r/| — > +00 and for Im(Y) = —71. 

Lemma 7.5 For all positive constants M, eq such that M > eq: 



Y\ + C 



□ 



hm l 3(t,^ 



ni(X - 1) 7/ m (y) 



_ 2i Y 

e * 



£ ln(2 V*V*v) 



-71 



2iy 
A - 1 



dY, 



uniformly for Im(r]/t 2 ) in compact subsets of (3/2, (3 + A)/2) and Eq < < M. 



Proof. We split the integral in (7.8) as follows: 

V2 



G t 



1L 

t 2 



iri (A 
V2 



1) / 



ImY=~Ti,\Re{Y)\< lr £-^ V(Y + T]/t 2 ) 



Virilt 2 ) _2 iy ( 2iY 
K " ' t wf 



■ t* I lnt| 



A- 1 



v^ri(A - 1) J ImY 



V(ri/t 2 



-^^<\MY)\< 



*o V(Y + v /t 2 ) 



r 



r „ / 2*y 



+ 



/Tri(A-l) 
Ji + J 2 + J 3 . 



Hv/t 2 ) 

\Re(Y)\> S -0 ViY + rj/t 2 ) 



t «r 



2iY 
A-l 



A- 1 



dy 



dy 



J\ = / [ — — - — t x ~ 1 — e A 

^i(X-l) J ImY =- lu \ Re (Y)\< i— r V^l^ + ^A 2 ) 

1 v ' 1 i z I m t\ 



A-ly v/wi (A - 1) i/ TO y=_ 71i 
_ ^_ / e -^m( 2 ^) r ( ^y A dy 

7Ti (A - 1) J/ m y = - 71i | jRe (y)|> V^- 1 / 
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Using that: 



e * 



i ln(^) 



iyl y l 



and 



2iY 
A-l 



<Ce ( A -!> 



-tx^tI^I 



the last term in (7.12) is estimated as: 



7d(A-l) J Im Y=- yi ,\Re(Y)\> 1 



_2.y 
e A 



£ ln(2 vW^) r 



2iy 
A-l 



dY 



< Ce t2 l ln 'l 



for some positive constant a. In the first term of (7.12), using (7.10) in Lemma 7.3 we are led to 
estimate: 



h(t) 



e (A-i) |r| |dy 



Im Y=— 7i 



which tends to zero as t — > + by (7.9). 

We now consider the integral J3 . This term is estimated using the estimates of the function V proved 
in Proposition 11.2 and Stirling's formula: 

I ^3 1 < Ctff «e i^U-i > <C s t 
for some positive constant a choosing 5 sufficiently small. 



The integral J2 is estimated using (7.11) in Lemma 7.3 

Vh\<c f 

J ImY=- 11>1T f^<\Re(Y)\<^ 

for some positive constant a. 
Lemma 7.6 

V2 f e -^ln( 2v ^^) r 



4(A-i)l y l e -A^rl y l|rfy| < Ce 



□ 



iri(\ - 1) i/ m (r)=_ 7l 



2iY 
A-l 



Proof. We first perform the change of variable: 2iY/(A — 1) = s: 



-V2 



ni (A - 1) Ji m (Y) 



-|i^ln(2v^v^) r 



-71 



2iY 
A-l 



dY 



1 



'2 it 



Re(s)=7i 



Then we deform the integration contour and reduce the integral to the sum of the residues of the 
function e~ s ln ( 2 V^m) a i th e poles s = —n of the Gamma function with residues - ^ . □ 
Proof of Proposition 7.2 We split the integral in (7.7) 



t 2 G(t, X ) =h + h + h 



where 



h 

h 
h 



1 



'2-k 
1 



Im(ri)=f3 t 2 , \ri\<eo 



e tr >XG[t,±)dTi 



e^G[t,^)d V 

2-7T Jlm(ri)=f3ot 2 ,Eo<\r)\<M V * ' 



1 



27T Jlm(ri)=l3 t 2 , \rj\>M 
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Using Lemma 6.1 

+h\<C f e^dr, < C (e + e~ aM ) 

Jlm( V )=/3 ^,\ri\<e ,\r,\>M 

The integral I2 is estimated using the previous Lemmas: 



{Git 



1 

'2-K J Im(ri)=l3ot 2 ,eo<\v\<M 

I 

J Im(r))=f3 t 2 ,e <\v\<M 



V 
t 2 



^IX e -2y^vl drl = J , 1 + J 2 2 



The first term is bounded as 

1^2,1 1 <C 



/m(r;)=/3ot 2 ,£o<l»?l<M 



Git 



t 2 



'2-KB 



and tends to zero by Lemma 7.5 and Lemma 7.6 And, 



lim I 2 a 



e i VXe -2V^n dr] _ 



Im(r))=0,e <\v\<M 



Now taking the limit as Sq —> 0, M — ► +00 and using (7.13) : 



lim t 2 G(t, X ) 



Im(ri)=0 



(7.13) 



(7.14) 



The integral in the right hand side of (7.14) may be calculated explicitly using contour deformation. 



For x < the contour is sent to the region where Im(r)) — > —00 and the integral gives zero. When 
X > the deformation is made to the upper half plane in such a way that it avoids the cut along 
the half line rj G i M + and the integral is reduced to : 

7T _ 3L 



1 J sin (2v-A) r-\\!,\ 



3/2 



e x . 



X 



□ 



7.1.2 Estimates of G{t,X) for t 2 < X < 1. 

We obtain an estimate in self similar form in the domain where t 2 < X < C that is required to show 
that G(t, X) converges to the Dirac measure as t — > 0. 

Proposition 7.7 For any 5 £ (0, 1/2), there exists a positive constant Cg such that 

4-2<5 



\G(t,X)\<C 



t 1 



x\l- 5 



for t 2 < \X\ < 1. 



Proof. We integrate by parts twice in formula (7.1) and obtain: 

d 2 



G(t,X) 



1 1 

l^x 2 



Im(0=Po 



Using that \e iX t - l| < C s \X\ 1 P +s \S\ 1 / 2+6 for 5 £ [0, 1/2], as well as Q we deduce: 

C<5 



|G(t,X| < 



X 2 



im(£)=A> 



(1 + lel 3 / 2 ; 



< 



Cat 



i^|i/2+^|i/2+* e -*vi€i|dei 
1*1 



Im(i)=f} 



1-5 



and the result follows. 



(7.15) 



□ 
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7.2 Behaviour of G(t, X) for t > 1. 

The behaviour of the function G as t — ► +00 has a self similar structure. This is seen by writing the 
function G(t,X), given in (4.3) and (7.1), in terms of the variable 

2 



(7.16) 



G(t,X) 



tt(A - 1) y/ m (0=A) 



Im(y)=/3 1 V(V) 



V(£) 



2/. 



A-l 



(£-1/) 



(7.17) 



with 3/2 < A) < A < (3 + A)/2 . 

Moving the integration contour of y downward in the expression of G(t, X), the first singularity to 
be met in the integrand of of (7.17) is y = i which is a zero of V{y) (cf. Proposition 11.3). That 
gives: 



G(t,X) = tA=r*i(0) + 



7r(A - 1) Jl m (O=l3 
*l(0) 



d£ e 



i£9 



Im(y)=/3 2 V \V) 

2i 



2i 



A-l 



(e-y) 



* V'(i) 



A-l 



(7.18) 
(7.19) 



A-l 

where now /3 2 6 ((3 - A)/2, 1). 

Notice that V'(i) 7^ by the following reason. If we differentiate the formula (4.5) we obtain: 

V '(») = -V *)*(-2~ *) " V (^~ 0* (-3- *) 
At £ = « the function V is analytic and has a zero. Therefore, 



By Proposition 11.3, V(^2^ i) ^ 0. Moreover, by Proposition 



10.17 



~ — 2 7T*(C — ^Y^ 1 *) as |i?e(C)| ^ +00 
and therefore i) = —2iri. We then obtain the expression for ^1: 



1 



2i 



A-l 



n(X-l)iV(^i) Jim(0=p 
We study now the function ^f(8) and give its behaviour as 6 — > ±00. 
Proposition 7.8 T/ie following estimates hold: 

M/ 1 (^) = C 2 e-^ + o(e-( A+1 - £ ) 
/or some positive constant e arbitrarily small and where 



(7.20) 



(7..S 



as 



-00, 
+00, 



(7.21) 
(7.22) 



C1 



2i V((l + |)i) 

ana 02 



(A-l) v(^i) 



2iri V( 



A+l 
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Proof of Proposition 7.8 . We use again contour deformation. In order to obtain the behaviour as 
— > — oo we deform the contour integration in downward. The first singularity of the integrand 



that we meet is £ = 3i/2 which is a pole of V(£). Using (4.5) and (10.17) we obtain: 



Therefore 



lies [ V, £ 
Ties ( V,£ 



3i 

(3 + A)i 



-tV 
V(2>) 

47TZ 



1 + - I / 

2 



(7.23) 
(7.24) 



2J V((l + A/2)i) 
(A-l)V((A + l)/2)i) 
1 

+ \— i 



e 2 C 



7r(A-l)iV(^i) Jlm(0=p 3 



d£e* V(f)r 



2i 



A - 1 



(f 



(7.25) 



r / A+l 



A-l J V(2i) 3 + \n 

e 2 



2vri V(^±i») 



+ 



tt(A- l)iV(^i) Jim(t)=ik 



2/ 



A-l 



(e-i) 



(7.26) 



where /3 3 G ((4 - A)/2,3/2), /? 4 G ((3 + A)/2, 1 + A). We have derived fl7^25| ), ( |7726] ) deforming the 
contour of i ntegr ation upward and downward respectively. 
Proposition 



11.2 



ensures that the function V(£)T ^ — ^r[(£ — ^ s integrable and then, for Re6 < 0: 



/m(f)=Ai 



2? 



A-l 



(€-0 



< Ce-P 3Re ( 0) . 



while for > 0: 



2/ 



A-l 



Proposition (7.8) follows from (7.25), (7.26), (7.27) and (7.28). 



(7.27) 



(7.28) 



□ 



7.2.1 Estimate of the remainder term in formula (7.18) 

Let us call that term: 



Gi(t, 



7r(A - 1) Jl m (£)=j3 Jlm{y)=l3 2 V{y) 



2i 



A-l 



(C-y) 



(7.29) 



We have then the following Lemma. 
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ImZ 



ReZ 



ImZ = 



Z c 



Figure 4: The curve T>\ 



Lemma 7.9 For any positive constant C\ there exist positive constants A, 5q and C 2 such that, for 
allt > 1 : 



\Gi{t,6)\e^ <C 2 t^~ &0 for all 6 <® 
\Gi{t,0)\e^ e <C 2 t^ &0 for all 9 >0. 

Proof of Lemma 7.9\ The function G\ is may be written as: 
Gi(t,0) = 



7r(A - J- J Jlm(Q=/3o 
Im(v)=Pa V (y> V A - 1 



(C-y) 



The function i?i is estimated in the same way as the function G(t, £) in Lemma 
first perform the change of variables: y = £ + and obtain: 



6.3 



(7.30) 
(7.31) 



(7.32) 
(7.33) 
To this end we 



m,o = VW\ [ 

J I 



dy 



v(0 



/m ( z)=tt v(^iei^ + e) 



a«(vieLg4j) 

-t x=i r 



2i 



A - 1 



lei ^ (7-34) 



Then we deform the contour of integration in (7.34) to the new contour T>\. Since along this new 
contour Im(Z) = , we have Im{ y/\£\Z + £,) < f3 2 and then 

2iy 2 ft, 2 _r 



We are then left with the term 



'Pi 



v(0 



t r - — V C ^ 



2i 



\dy\ 
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that may be estimated following the same arguments as in the proof of Lemma 6.3 We obtain then 
the bound on Hi(t,£): 



(7.35) 



Using (7.1) we deduce that, for any positive constant R there exists a constant Cr > such that, 



for all 101 < R: 



\Gi(t,0)\ < C R t&- S 



for some positive constant Cr. This shows ( 7.30[ ) and (7.31) for 9 bounded. 

In order to prove (7.30) and (7.31) we deform the contour of the integral with respect to £ in 



singularity of G\(t, £) is located at £ = (3 + A) i/2 (see Proposition 11.3). Therefore: 



formula (|7.29|). To obtain the estimate (|7.31|) as 9 — > +oo we defor m the contour upward. The first 

(7.36) 



3+A, 



Gi(i,0) = b 2 (t)e-^° + Q 1 (t, 



Qi(t,9) 



- 1) Jlm(Z)=f3 5 Jlm{y)=Pi ^KV) V X ~ l 



it-y) 



(7.37) 



where (3 5 G ((3 + A)/2, (4 + A)/2), and 

i /,\ V(2i) / , / 2i /3 + A. 

Using that again (3\ G ((3 — A)/2, 1) and so in particular (3\ < 1, we deduce: 

\b 2 (t)\ < Ct^~ 5 °, for t > 1. 
We are then left with the term Q\(t, 9) that we write: 



Qi(t,e) 

H 2 



i 



V(£) at* 



2i 



(7.38) 

(7.39) 
(7.40) 



vr(A - 1) J Irn{y)=f3l V{y) V A - 1 

The function H 2 is estimated in the same way as the function G(t, £) in Lemma 6.3 We change of 



variables: y = £ + y |£| Z and obtain: 



H 2 (t,0 = VW\ 



dy 



V(0 



2i(-/|?Tz+o / 2i /— r 



(7.41) 



Then we deform the integration contour in (7.41 ) to T> 2 Since along this new contour Im(Z) = J^* ; 



we have Im(\J\(\Z + £) < (5\ and then 
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ImZ 



ReZ 



ImZ = 



z c 



Figure 5: The curve T>2 



We are then left with the term 



©1 



v(0 



2i(yjg|g+g) 

-f a=i r 



2i 



A- 1 



that may be estimated following the same arguments as in the proof of Lemma 6.3 We obtain then 
the bound on H2(t,£): 



Using now (7.39) we deduce that: 

\Qx(t,0)\ < Cf&-- S e-M. 



(7.42) 



(7.43) 



Combining this with (7.38), estimate (7.31) follows. 



The estimate (|7.30|) for 6 — > — oo is obtained in a very similar way. We deform downward the 

of 
□ 



contour of the integral (7.33) and we continue the proof as for (7.31). This concludes the proof of 
Lemma 17.91 



7.3 Behaviour as < t < 1, \X\ -> +oo. 

For the small values of time the solution is described in the (t, X) variables as follows. 
Lemma 7.10 For < t < 1 the following estimates hold: 

' e- S 2 X t + o(e~^) x t) as X 



G(t,X) 



-oo 



where: 



Cl e-^ x t + o(e~^ +£ ) x t) as X^+oo. 
V(2i) 



(7.44) 



Ci 



4vrV((l + f) i) 
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Proof of Lemma 7.10. Using (7.1) we obtain: 

e* x G(t,£)d£ 



G{t,X) 



1 

fair Jim{Z)=p 
i 



7r(A - 1) Jlm(0-- 



[ e* x dt [ 



dy- 



v(0 



m( y)=- 7l V(Z + Y) 



2i v 



2iY 
A-l 



where /3o £ (§, ^2 ) and 71 small. Lemma 
for large |£| and this yields convergence of t 
for f e (| - So, ^ + <5 ) with 5 > 0. 



6.1 



implies that G(t,£) is exponentially bounded on |£| 
re integral. Moreover, such exponential decay holds also 



We can now deform the contour integration on £ crossing the poles of G (t, £) that are due to 



the poles of V(£). The closest poles are at £ 



2 ■ 



£ = ^r^« respectively. We deform the contour 



upwards if X > and downwards if X < 0. We then obtain using (7.23), (7.24): 
G{t,X) = 



-2vr 



vr(A-l) 



V 



1 + 



A 



/ I e~^ x 



dy- 



t A-i- 



/ m (y)=-7i V(fi + y) 



2iY 
A- 1 



+ 



+ 



7T(A - 1) Ji m (£)=fa 
Jl + J 2 



d( 



dy 



v(0 



7m(Y)=- 7l 



G(i,X) 



7r(A - 1) Jl m (£)=fo 
= Jl + J2 



e 2 



V(2i) t 
2vr(A- 1) 



_2j_ V 



Y 



— r 

/m(y)=- 7l V(^i + Y) 

v(0 



Jm(y)= _ 7 / y V(£ + F) 



2j y 



2iY 
A-l 



2iY~ 
A-l 
2iY 



(7.45) 



A-l 



(7.46) 



where /3i 



5 , /3 2 = 3±A + j w ith 5 > small. 



The terms J2 in (7.45), (7.46) can be estimated easily. Indeed, they can be written in the form: 

1 



G(t,X)~- 
Integrating by parts we obtain: 
G(t,X) = 



27T Jlm(Z)=/3 e 



e* x G(t,Z)d£ , £ = 1,2 



^X 2 



die 



Im(0=l3 e 



G(t,0 



We can now estimate 



(G(t,0 



using Lemma 



6.1 



It then follows that: 



U2I < 



C_ 
X? 



1 




>Im{£)=P t 





-ay/\i\t 



3/2 



Ce^ x 
X 2 



-WW 



Im(0=Pt 1 + |£| 



,3/2 



1*1 



In order to estimate the integral we split it as follows: 



I, 



(i + iei 



-ay/\£\t 



3/2 



[ [-] m + / 



< Ct + Ct 2 



[...) \dt\ 
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Then: 



I J2I < Cte^ eX for \X\ > 1 



where I = 1 for X < and 1 = 2 for X > 0. In order to compute the terms J\ in (7.45), (7.46) 
we deform the contour on Y upwards. We cross the first pole of the function T ( l^j ) for Y = 0. 



However, this point is not a pole of the integrand, because the functions V(fi+Y), V(^i + Y) have 
also a pole at Y = 0. Therefore the first pole of the integrand that is found is the one at Y = ^ipi. 
Notice that V(|i + 1") does not have a zero before, since 1 + = l + f < 1 + f ■ On the other 
hand ^ + ^ = 1 + A while the first zero of V (1?) is at 77 = 2 + | and 1 + A < 2 + |. Then, after 
deforming the integral contour as indicated, the terms J\ can be written as: 

ay 2i , // AN N 3x f , „( 2iY 



and 



A-l VV 2 7 / " Jlm(Y)= l2 V{\i + Y) \X-1 
e"i x (t+J 3 ) 



_3+Av V(2i)t 
Ji = e 2 



47rV((l + f) i) 

iV(2i) _3±a x r v( 2iY 

~2vr(A-l) e 2 A m(y)=72 y V(^i + y) VA^T 

where 72 > (A — 1) /2 is such that 72 — (A — l)/2 is small. In both cases there exist positive constants 
5 and C$ such that for all t > 1: 

|J 3 | < C s t 1+5 . 

\Ji\ < Cte PeX for \X\ > 1 

where £ = 1 for X < and £ = 2 for X > 0. 
It then follows that: 

|Jl| + IJ2I < Cte^ eX for \X\ > 1 



where, as before, 1=1 for X < and ^ = 2 for X > 0, and Lemma 7.10 follows. 



8 The initial value problem. 



Using the fundamental solution g obtained in Theorem 2.5 we can obtain a solution of the initial 
value problem 



dh 
h(0,x) 



L[h) 
h (x) 



1.2) 
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with L[- \ defined in (2.2 ). Assuming that there are not difficulties with the integrals written below, 
we would expect, due to the linearity of the problem (8.1 ), (8.2 ) the following representation formula 

dy 



for their solutions (cf. Theorem 2.5): 

h (t, x) 



ho (y)g ( ty x ^ , -,1 

y J y 



(8.3) 



We first precise sufficient conditions on ho that allow to define h(t,x) in (8.3). 
Theorem 8.1 Suppose that the function ho £ C satisfies 



Jo 



ho{y)\y X dy+ I \h (y)\ dy < oo 



•4) 



Then the function h(t,x) defined for t > 0, x > by means of (8.3) solves the initial value 
problem \8. i| ), \8.3j) . 

The proof of this Theorem reduces to a detailed analysis of the conditions on ho yielding integrability 
of the right hand side of (8.3). 

Under more stringent assumptions on ho it is possible to use Theorem 2.5 to derive more detailed 
information on the asymptotics of the solutions of (8.1 ), (8.2 ) for x — * and x — > oo. The meaning 
of this asymptotics will be explained in the next Section. 

Theorem 8.2 Suppose that 



\h (x)\ < Cx~2 +£ , 0<x<l , e>0 
\h {x)\ < Cx-( 1+£ ) , x > 1 , e > 



Then the function h(t,x) given in (8.3) satisfies for any t > 



h(t,x) -A- (t)x"i < (t)x~^ +e for < x < 1 
h {t, x) - A+ (t) x~^r < B + (t)x~^r~ e for x > 1 



J.5) 
1.6) 



(8.7) 



for suitable functions A_ (t) , A + {t) , S_ (t) , B + (t) . 
Detailed proofs of these two results will be given in [11 . 



9 Particle fluxes for singular solutions of the coagulation equation. 

9.1 Computing particle fluxes for the nonlinear equation. 



It is well known that smooth solutions of the coagulation equation (1.1), (1.2) preserve the value of 
the quantity J Q °° xf (x, t) dx. This could be expected given that the equation can be interpreted as 
a description of the coalescence of particles with sizes x, y to form a particle of size (x + y) in the 
time interval [t, t + dt] taking place with a probability K (x, y) dt for each particle pair (x, y) . 



Notice that, in spite of this conservation property, we cannot expect to be able to rewrite (1.1 ), 



(1.2) as a conservation law with the form 
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with j depending only on local properties of / at the point (t, x) due to the fact that the particles 
can change their size an amount of order one in an infinitesimal interval of time [t, t + dt] . However, 
we can expect to have equivalent ways of finding formulas with the form: 



R2 
Ri 



jil xf(x,t)dx)=J+ uR2 -J RiR2 (9. 



for arbitrary values < R\ < R2, where J R R ^ denotes the number of monomers that enter in the 
interval x G [R\ , R2] coming from particles with sizes in the region IR + \ [R\ , R2] for unit of time and 
J Ri R2 is the number of monomers that leave the set x £ [R\, R2] for unit of time. 

A careful counting of the number of particles entering and leaving the interval [R\, R2] yields 

j£ lRa = \l f(y)f(z)K(y,z)(y + z)dydz + 

+ [ f(y)f(z)K(y,z) ydydz (9.10) 

JU 2;R 1 ,R 2 

D t; Rl ,R 2 = z):y<R 1 ,z<R 1 ,R 1 <(y + z)< R 2 } 
D i Rl ,R 2 = i(v> z ) ■■ V < R h Ri < z < R2, Ri < (y + z) < R 2 } 



J 



+ 



D 
D 
D 



2;Ri,R 2 



3;Rx,R 2 



L 

{(y,z) 

{(V,z) 



1 

2. lri 



f(y)f(z)K(y,z) (y + z) dydz + 



/ (y) f (z) K (y, z) ydydz + ! f (y) f (z) K (y, z) ydydz (9. 

2;R 1 ,R 2 D Z\R 1 ,R 2 

■ Ri < y < R2, Ri < z < R 2 , R 2 <(y + z)} 
:Ri<y<R2,R2< z, R 2 <(y + z)} 
:Ri<y<R2, z<Ri, R 2 <(y + z)} 



11) 



There are two particular cases of (9.9) that will be particularly relevant for our purposes, namely 
the cases 

R 1 = , R 2 = R G (0, 00) (9.12) 



and 



Ri = Re (0,oo) 



Ri 



00 



In the case (9.12), (9.9) reduces to: 



d_ 
dt 



R 



xf (x, t) dx 



-Jr [/] 



).13) 



).14) 



J R if] 
D~ (R) 



D l:R 1 ,R 2 



f (y) f (z) K (y, z) (y + z) dydz + ! f (y) f (z) K (y, z) ydydz 



{{y,z):y<R, z < R, R < (y + z)} 
{(y,z):y<R, R < z, R < (y + z)} 



(9.15) 
(9.16) 
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In the case (9.13), equation (9.9) becomes: 



dt 



xf(x,t)dx) =Jt[f] 



R 



).17) 



41/] 



/ (y) f (z) K (y, z) (y + z) dydz + f (y) f (z) K (y, z) ydydz 

'Dt(R) JD+{R) 

D+{R) = {(y,z):y<R, z<R, R < {y + z)} 
D+{R) = {(y,z):y<R,R<z,R<(y + z)} 



Using the formu 
the solutions of (1.1 



a ( 9.14 1 we can easily prove that for the kernel K{y,z) = (yz) 2 the meaning of 

1 is that there is a flux of particles 



1.2) having the asymptotics f (x) 



3+A 



escaping to infinity. Indeed, it is readily seen that for this kernel: 



lim J' [/] = 2ttA 2 

R^oo 



Notice that the function / (x) 



3+A can be thought as a stationary solution of (1.1), (1.2) charac- 

X 2 

terized by a constant flux of particles propagating from the smaller to the larger values of x. 



9.2 Particle fluxes for the linearized equation. 



Since equation (2.1), (2.2) has been obtained linearizing (1.1), (1.2) we can derive formulas for the 



particle fluxes associated to (2.1), (2.2) linearizing (9.9 ). We recall that (2.1), (2.2) has been obtained 
using that 

.18) 



3+A 



f(x,t)=x 2 +g(x,t) 



in (1.1), (1.2) and keeping just linear terms on g. Plugging (9.18) into (9.12) and keeping just the 



linear term on g we obtain: 



dt 



xg (t, x) dx 



-J 



R,lin 



(9.19) 



Where J R,lin is: 



J 



R,lin 



Ds(R) 



L 



^Tfx(yz)Hv + z)dydz + 



+ 



D-(R) y~2 

9(z) 



3+A 

D-(R) y 2 

A 
Z2 



— 

(yz) 2 zdydz + 



^l(yz)2(y + z)dydz 
D-(R) z~ 

^^(yz)^ zdydz 



D-(R) z 2 



Z 2 



-1 



a(y + z)g{z)dydz + / -^g(z)dydz + / -^g(y)dydz 
Da(R) y* Jd^(R) y2 JD^(R) Z 2 



y 2 



h(t,R) + I 2 (t,R) + I 3 (t,R) 
{0<z<R, < y < R, R < (y + z)} . 



and D 1 {R) ,D 2 (R) defined by (9.15) and (9.16) respectively. Integrating (9.19) we obtain: 



R pt r-R 

xg(0,x)dx= / JRu n ( s )ds+ / xg(t,x)dx 



(9.20) 
(9.21) 



.22) 
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The solution g(t,x) obtained in Theorem |2.5| satisfies 

g(t, x) ~ a{t) x~ (3+A)/2 as x -> +00. 



).23) 



Taking the limit i? — > 00 in (9.22) and using the fact that the terms I\ (t,R) , 1% (t,R) tend to 



zero as R — > 00 and that lirnR—xx, I2 (t, R) = a (t) , it follows that: 



00 rt poo 

xg(0,x)dx = / a(s) ds + / xg(t,x)dx. 
Jo Jo 



).24) 



The self-similar asymptotics (2.10)-(2.13) implies that lim^oo xg (t, x) dx = 0. Then: 

00 POO 

xg(0,x)dx= / a(s)ds. 
Jo 



).25) 



The left hand side of (9.25) is the initial total mass of the perturbation. The right hand side of 



(9.25) is the total amount of particles contained in clusters of infinite size. Equation (9.25) means 



that all the excess of particles initially introduced in the system move as t — ► +00 to an infinitely 
large cluster. 



10 Appendix I. 

We take now the Mellin transform on both hands of the equation. We recall that the Mellin transform 
of a function g(y) is defined as: 



M(h)(s) 



y s 1 h(y)dy. 



(10.1) 



Taking the Mellin transform of the right hand side of (2.1), (2.2) we obtain after straightforward 
calculations: 



0_ 

ds 



M(g)(s) 
M(s) 



51/2-s 



l)-3/2_0-3M d9 + 

2- s 



+ / (1 - #r 3/2 (6 s - 1 -l)dO+- 2y/2 

Jl/2 s 

h(s) + h(s) + h(s) + h- 



(10.2) 



(10.3) 
(10.4) 



Remark 10.1 // the function g satisfies the estimates (2.5), (2.6) for some r > 1 and p < 2 such 



that p < r we will have that its Mellin transform is well defined in the strip Res £ (X/2 + p, X/2 + r). 
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10.1 The function M(s). 

We consider in this Section the auxiliary function M obtained by taking the Mellin transform of the 
equation (2.1), (2.2) and first rewrite it in terms of Gamma functions. 



Proposition 10.2 The function M(s) defined in (10.3) can be written as 

M(.S) : 



2y^r(s) 
r(--i/2) ■ 



(10.5) 



Proof. We can write the term I\ of (10.4) as 

h{s) = 



l/2- S 0-3/2 J f 1 



-3/2 



1 \d6 



Using Binomial Theorem to expand (l — i) 3//2 and integrating each term of the resulting series we 
obtain: 



71=1 



-3/2 
n 



■iy 



2~n 

s + n 



Then, 



-3/2 
n 



-iy 



2 -n 



s + n 



/ 1 + / 3 = 2- s ^ 

n=0 

Integrating by parts we obtain: 

I 2 (s) = 2V2-2V2(l/2) s ~ 1 -2(s-l) / X (1-^)~ 1/2 ^ 2 ^. 

Jl/2 



(10.6) 



(10.7) 



In order to compute the last term in (|10.7|) we write: 
»i 

1/2 



pi r l/2 

(i - e)- l/2 e s - 2 de = / (i - e)- l i 2 e s - 2 do - (i - 9)- l i 2 e s - 2 de, (10.8) 

io io 
Expanding (1 — #) -1 / 2 using the Binomial Theorem and integrating each term of the resulting series: 



1/2 \ (-i) m _; 



Moreover, 
2 



1 T \ n + i 



n=0 

where we have used that 



-1 



in+l 



n + s 
-1/2 



n + 1 



2(1 ~ 1/2)- 
s - 1 

-3/2 



-1/2 



2( S -1) 



n=0 



1)^U+1 



3/2 \ (-l) n 



n + s 



(10.9) 



2~ n (10.10) 



(10.9) and (10.10) we deduce 



1/2 



1 - 9)- 1/2 s - 2 d6 = 2" 



2(«+ 1 ) y n 
, f 2^2 



as well as the Binomial Theorem. Combining 



s-1 2(s-l) 



1) ^ ( 71+1 J 
7 n=0 v 7 



3/2 \ (-ir 



n + s 



Therefore, using ( |10.3[ ), fll0.6| ), (jlO/Tj), ( |10.8| ) and ([10T1J, 

M{s) = -2{s - 1) / (1 - 6)- l l 2 e s - 2 d6 



2 v / 2T(s) 

'r(«-i/2)" 



(10.11) 

(10.12) 
□ 
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10.2 Zeros and Poles of M(s) 

The function M (s) is meromorphic on the whole complex plane. Its zeros are given by the poles of 
the function T(s — 1/2): 



s z(n) = \~ n ' n = 0, 1,2, • • • 

The second zero, M(— 1/2) = corresponds to the solution x~^ 3+x ^ 2 . 
On the other hand, the poles of M{s) are those of the function T(s), i. e. 



(10.13) 



s p (n) = —n, n = 0,1, 2, 



10.3 Fourier variables. 



We make now a change of variables in order to have functions defined in all of the real line M. To 
this end we define: 

M(g)(s)= g(y)y s - 1 ds= g{e x )e~^dX = V2^G{£) 

JO J-oo 



i.e., M(g)(s) = V2ttG(£,) with s = then 

M(g) (s + ^ 



'2ttG U + ^-i 



If we re-write the equation (10.2) in these new variables we obtain 

It is more convenient to use instead of M the function defined as: 

$(0 = M ( if + 1 + A 



Using (10.5), we obtain <]?: 



HO 



2y¥r(ie+i + t: 



The equation (10.15) reads then as equation (2.20). The set of zeros of $ are: 



The set of poles is 



Notice that by definition, 



f z (n) = i (n+^ti ] „-o. |.2.--- 



n = 0,1,2, 



i p (n) = i + 1 + 



M 



1 



I I 1 • A / i -o. 



(10.14) 



(10.15) 



[10.16) 



(10.17) 



[10.18) 



[10.19) 
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10.4 Behaviour of $ at infinity. 



Theorem 10.3 Let us define = (s) = sgn {Im (s)) . For any fixed L > 0, the following asymp- 
totic formulas hold: 



M(s) 
M(s) 



-\ > i 1 - f +() ' 

8s 



as |/m (s)| — > oo 



n 2m (1 + i9) V-isG - i — + 



3/2 



(10.20) 

as \Im (s)| — ► oo 

(10.21) 



uniformly in sets where arg (s) G (— ir + £q, vr — eg) / or aK 2/ £ o > 0. 



Proof Formula (10.20) is a consequence of (10.5) as well as the asymptotic formula: 



12z \z 

that is uniformly valid in sets arg (z) G (— ir + eq, it — eq) for any eq > 0. 
Then: 

M(s) 



as \z\ — > oo 



(10.22) 



as s — > oo, 



uniformly in sets arg (s) G (— 7r + £o ; tt — eo) an< f arg (s — 1/2) G (— ir + Eo, it — Eq) for any eo > 0. 
Notice that: 



i 



as s — > oo, 



uniformly in the same sets as above, whence: 
M{s) = 



2 v '~.s| l-|- + ' 

8s 



as s — ► oo, 



uniformly in sets arg (s) G (— ir + eo, t — £o) for any eo > 0. and (10.20) follows. Formula (10.21) is 
a consequence of (10.20). □ 



11 Appendix II: Some technical propositions. 

We prove now some auxiliary results used to prove the results of the paper. 
Lemma 11.1 Suppose that f and h are two analytic functions in the cone 

C{2e ) = {C G C; C = \C\e ie , G (-2e ,2e )} 

for some Eq > and real valued in M + . Let us also assume that 

\f{re i9 )\ + \h{re i9 )\ 



1+r 2 



-dr < +oo, for any 9 G (— 2eq,2eq) 



lim /(C) = 0i and lim /(C) = 2 



\f'(0\ = o 



as |C| -0, |CHoo, CeC-(2e ). 



(11.1) 
(11.2) 

(11.3) 
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Then, the function 



1 

2tH 



(h(s) + if(s)) 



C 8 + 1 



ds 



is analytic in the domain: 

D(e Q ) = {C 6 C; C = \C\e W , € (-e ,2vr + e )} . 

Moreover, 

9i 



no 



2tt 



lnC + ^(r?)+o(ln|C|), as ( -» 0, 



02 

"2tt 



lnC + ^(r?)+o(ln|C|), as |C| -> +oo, (6%). 



where the function H(£) is a real valued function defined by 

1 1 



#(C) 



1 

"2^ 



C s + 1 



(is. 



(11.4) 
(11.5) 

(11.6) 
(11.7) 

(11.8) 



Proof of Lemma 



11.1 



Using the Lemma C.2 of [9] with the function / we obtain (11.6) and 
(11.8). On the other hand, the condition ( |11.1[ ) ensures that the function F is well defined and of 
course real valued. □ 



Proposition 11.2 Let V(£) defined by (5.6) and (7.1). Then, for any 5 > arbitrarily small and 
all M > arbitrarily large, there exist two positive constants C\ t $ m,C2,5 m such that 



y 2 . 



h5)]£| 



< |V(0I <C 2AM e 



-(-- 

V 2 . 



(11.9) 



uniformly for Im(£) in compact sets of (3/2, (3 + A)/2) as well as for all £ such that \Re(£)\ > 1, 
|Jm(f)| < M. 



Proof of Proposition 11.2 Given £ such that Im(£) G (3/2,(3 + A)/2) we can represent the 
function V by (5.6) and (7.1) for (3\ such that (3\ — (A — l)/2 < £ < j3\. In order to simplify some of 
the calculations we use the following change of variables. 



KC) = v(£) 
v(C) = *(0 



(11.10) 

(ii.il) 

(11.12) 



The function ln(— tp(s)) may be written as 

In (-¥>(*)) =ln(|^(s)|) + iarg(-^(s)). (11.13) 

The functions ln(|</?(s)|) and arg(<£>(s)) satisfy the hypothesis required to h and / respectively in 
Lemma 11.2 In particular, by Proposition 3.1 and the fact that Re{— 3>(y)) > for all y such that 
Im(y) 6 ((2 + A)/2, (3 + A)/2) we may normalize the argument of the function In (— y?(s)) such that: 



limaxg (-<?(£)) 



7T 

'4' 



7T 



lim arg(-y>(0) 
4^00 4 



(11.14) 
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Applying Lemma 11.1 it follows that: 

1 



1 

2ni 

1 

2vrf 



In ("¥>(*)) 
In ("¥>(*)) 



1 



S-C 8+1 
1 1 



S-C 8+1 



ds 



ds 



1 



ln(-C) + i-H"(r?) + o(ln|C|) 



(11.15) 



as C -» 0, C G O(eo) 
1 



ln(-C) +iH(rj) + o(ln|C|) (11.16) 
as (" — » oo, C G -D(ffo)- 



The two estimates in (11.9) follow, for Im{£) in compact sets of (3/2, (3 + A)/2), by taking expo- 
nentials in both sides of (11.15) and ( 11.16| ) and inverting the change of variables (11.10 )-( 11.12 ). 
In order to prove the estimate for £ in the region |i?e(£)| > 1 and |/m(£)| < M, we extend analyt- 
ically the function V(£) to such regions using (4.5) as well as the fact that, by Proposition 3.1 we 
have for some positive constants C\ and 0%: 



cM l2 < l*(OI <«| 1/2 

for \Re{£)\ > 1 and \Im(£)\ < M. 

Proposition 11.3 The zeros of the function V are: 

1 + - + n + k- — 1 | i, n = l,2,---, k = 0, 1, 
2 2/ i)j 

1 + A ,A-1\ 

k — - — | i, k = 1,2, 



The poles of the function V are: 

1 + A 



+ n + fc^; — ) i, n = l,2,---, = 0, 1, 



1 A ,A-r 

H k i i, 

2 2 



A; = 1,2, 



□ 

(11.17) 
(11.18) 

(11.19) 
(11.20) 



Proof. The proof follows from (4.5) and the distribution of zeros and poles of $ given in (3.1) and 
(3.2). Consider first any ij 6 C such that / 771,(77) > (3q + (A — l)/2 and let k be the first natural 

Im(rj) - k^—^ <G [Po,Po + 



integer such that 



By (Of, 

Vfa) = (-ir 

as far as none of the points 77 



v(q-j^) 

$(77)^(77 - ^=li) • • • $(7/ - J^i) 



A=li 
2 



0, • • • j is a zero nor a pole f <£. Therefore, if 
.A-l 



Vj = 0,1,- • -k : r]-j——i {6(n),f p (ra), n = 0,1,- • •} 



then 



V(r?) = (-1)* 



$(77)^(77 - ^-M) • • • $(77 - kK^i) 



(11.21) 



;n.22) 



;n.23) 
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and therefore V(rj) is well defined and no t zero. If for some j E {0, 1, • • • k} and n E {0, 1, • • •} we 
have i] — k^^i = £ z (r&), then, by (11.26), the point r\ would be a pole of V. Moreover, that is the 



only way for a given r\ to be a pole of V. Therefore, the poles r\ of V such that lm{r\) > (3q + (A — 1)/2 
are: 



1 + A ,A-r 

n H h k — - — | %, 



n = 1,2, • --, A; = 0,1, 



[11.24) 



,A-1, 



f p (n), then, 

by (11.26), the point r/ would be a zero of V. And this would again be the only way for rj to be a 



If, on the other hand, for some j E {0, 1, • • • k} and n E {0, !,-••} we had rj — k'- ± ^ ± i 



zero of V. We deduce that the zeros rj of V such that Im(rj) > /3q + (A — l)/2 are: 

n = 1,2,---, A: = 0, 1, • • • 



A ,A-r 

1 + n H h A; I 

2 2 



(11.25) 



We have thus proved that, above the strip Im(rj) E (/3q 5 /So + (A — 1) /2), the zeros and poles of V are 



given by (11.17) and (11.19) respectively. 



A similar argument may be done when Im(rj) < Pq. If A; is the least natural integer such that 

Im{r)) + A;^^ E (A), Pa + 



we have by (4.5), 



vfa) = (-iy 



$(7?)$(ry + ^= i i)---$(r/ + i^= i «) 



(11.26) 



as far as none of the points rj + £^^i, £ = 0, • • • j is a zero nor a pole f $. Notice that by the 
definition of k and the choice of (3q, the only zero of $ that my be founded is £ z (0)) = (1 + X)i/2 
and the only pole is £ p (0) = (2 + \)i/2. The proof then follows as before. □ 



12 Appendix III: Stationary phase. 

We must estimate in Sections [6] and [7] several integral expressions of the form 



ImY=f3 -Im£ 



for a given function \& but different functions A. This is done using the stationary phase method. 
We collect in this Section some technical results about the function Y, t) and its critical points. 

12.1 The critical point. 



We first compute the critical points of the function Y, t) defined by (6.34) or equivalently, that 



of $(£,Z,t) defined in (6.10). Define: 



D(£, B) = \ Z E C; ImZ < 0, \ImZ\ < B 



ReZ+Qy^l 



sign(ReZ) = sign(Re^) > (12.1 



where Q = sign(Re(£)). 
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ImZ 




ReZ 



Figure 6: D(^,B) when Re(£) > and then Q = 1. 



12.1.1 The case < t < 1. 

We start with the following: 

Lemma 12.1 Consider the function F(£,Z) defined by means of 

2 



(A — 1) i Jlm^^fc 

for fa G (3/2,(3 + A)/2), Im(f) G (3/2, ft) and 

1 (A-H)) A_l 



ln(-cl>(n))e(7 ? -e,v / l^l^)^ 



< Im(Z) < 



(/?! - 7m(0) 



(12.2) 



(12.3) 



For any constant B > f/je function F(£, ■) can be extended analytically on the variable Z to the 



domain \Z\ < jp-, G D(£,B) if |£| > £o / or Co = Co(-S) sufficiently large. Moreover, there exists 
a positive constant C = C{B) such that 

21 



for \Z\ < 



ZeD{£,B) and |£| > 



< C Z 2 + 



(12.4) 



Proof. The function F is well defined in (12.3) since in that domain the variable Z is in the 
region where the function is analytic. The function F(£,-) given by (12.2) is then analytic in 
the strip \ImZ\ < ~^j|p f° r some 3o(I m O sufficiently small. We now claim that for any fixed 



constant B > 0, and Im£ G (A),/3o + (A — l)/2), if Re£ is sufficiently large (depending on B), this 
function F(£, •) may be extended analytically to the region.D(£, B). To prove this we derive new 
representation formulas of the function F performing suitable contour deformations in the variable 
7] in Q12.2D. 



Notice that the singularities of the integrand in (12.2) are contained in the set Rer] = 0, r\ = 
i + (A - 1)^/2 and rj = £ + s/W\ z + (A - 1)V 2 f ° r 1 £ z - Given Z in the region above, let be Z 
such that: 



i?eZo = ReZo, \ImZo\ < 



5 (Im£) 
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Imij 

















> > 






M\zf 


> 

Im rj = 






















M\Zo 





Figure 7: The curve C2 



and the integral curve Im rj = /?o + — £ lies between the two points £, + Zq and £ + Zo + 
(A — l)i/2. (Notice that this is possible since So(Im^) may be made as small as we need and 
Imrj > I ml;.) 

Consider the vertical segment of the complex plane connecting Zq and Zq: Zq = (1 — 9)Zq + 9Zq, 
9 £ [0, 1]. We then obtain an analytic extension of F(£, ■) varying 9 continuously from to one and 
deforming continuously the contour Imr] = (3q + — e in such a way that: 

• it always passes between £ + ZqsJ\£\ and £ + Zq + (A — l)z/2 



• we do not change the original integration contour for 

\MZo)y/W\ 



Re (p-£-Z y/W\) 



> 



The first condition ensures that the integration contour never crosses any of the singularities of 
the function ^1 — e~^ [ ^ , ^^ Z ^'^ . The second one ensures that it does not cross neither any 

of the singularities of ^1 — e^ 1 ^'^ 

Finally, since sign(ReZo) = sign(Re^), the new integration contour never crosses the line Rerj = 
where the singularities of In (— <&(t/)) are located. 

To estimate this integral we write 

^— I In (-$(,,)) 8fo - g, VW\Z)d V = 2 / &(r ] -tVW\Z)d V + 

~ l ) 1 JC-2 l A ~ L ) 1 JC 2 



(A 
+ (A 



L ln (m) @{v " e ' ^ z)dv =h+ h ' (12,5) 
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In (-<*>(£)) y/\Z\Z. 



(12.6) 



The first integral, I\ is computed explicitly : 

_ 2 1n(-$(Q) m 2i 

In order to compute the second integral we have to distinguish the cases i?e£ — ► — oo and i?e£ — > +00. 
Since both may be treated using similar arguments let us treat only the case i?e£ — > +00. In that 
case we decompose 1% as follows 

2 



h 



(A — 1) i JC 2 ,Rer l >0 ) \r)-Z\<^ 



■}dr] + 



+ 



(A-l)i 



/ [■•■]*; = / 2 ,i + J2.2 + k 

JC 2 ,Rerj<0 



C 2 , J Re 1? >0,|T7-e|> 1 | i 
2,3 



[•••] dv 



(12.7) 



In l2 t i we use now Proposition 3.1 and Taylor's expansion to obtain: 



$(77) _ V 7 ? 



using now 77 = £ + (77 — £) we can write: 



Therefore: 



In 



$(77) 

*(0 



/<&(7?)\ 



1 + 



v - £ 



+ 



as ,Re£ > 1, 1 77 - £| < 



l£l 



c) [ "LA ) + o 



iV asi?e£>l, |r7-£|<M. 



We now estimate the two following integrals for all Z £ -D(£, 5) and |Z| < -^5—. The first one: 



. , . ^) 6(77 - £, ^Z)d77 = I 



C 2 ,Rea>-Re£, \a\< 1 -^ 



< 



l£l JC 2 ,Rea>-Re£, W\<2\^\ 1 / 2 Z 
1 

+ 1 



< CZ 2 + C- 



a9(a,Vl£|^) 



r/cr 



*e(a,y/\£\Z) 



l£l Jc 2 ,Rea>-Re£, |<r|>2|5| 1 /2z 



(2(7 



-o[C|V3|Z| 

"lei 



for i?e£ > £ (5) 



the second one: 



L 



C 2 ,Rer,>0,\ V -t;\<^ 



¥ l£l 2 



e(7?-£,Vl£l^) 



drj 



2 /- 



1 



C 2 ,Reo->-fieg,|cr|<J|l 



(id + 



e(<W|£|z) 



(id 



C 2 ,i?e ->-Re^,| -|<2|^| 1 /2^ 

+ — / 

kl 2 ic 2 ,-Re CT >-ReC,|cr|>2|C| 1 /2z 



©(^viei^) 



ricr 



* ^ + - c \} z]+ w) + C ^r~ for fee > UB) 
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where in both cases C2 = C2 — £j £o(B) is a positive constant sufficiently large, depending on B and a 
and C are positive constants which may depend on B but are independent on Ret; and Z £ -B). 
From where we deduce that, for all i?e£ > ^q(B): 

e -a\Z\ 1/2 \Z\ 

- ^- (12.8, 

In the integral ^2,2 we use the fact that when \rj — £| > ^ and Re£ > £o(B), the function 
®(v ~ £> ^0 nas an exponential decay C7e _a ^' with C and a as above as well as the inequality 
I ln(-<J>(r/))| < C\ ln(\r] - £| + |£|)| < C|(ln(|r/ - £|)| + | ln(|£|)|) for 77 is large. For 77 of order one we 
use that ln(— $(r?)) is of order one to derive a similar estimate. Then 



\h,\^c[\z\ 2 + j^)+c- 



'2,2 



< 



3 -a|»7-fl 



(|ln(|r/-e|)| + |ln(|el)l)^ = 0fe-^ 1 



(12. 



'l*?-«l>f 

Finally, the estimate of I2 3 follows using the same cut off properties of the function O since Rer/ < 
and Re£ — ► +00 implies that Re{r\ — £) > C|£|. The final estimate of I2, by (12.7), is then, 

e -*\a\ 1/2 \z\ 



\h\<c\\z\ z + 



1 

w 



+ c- 



Using Proposition 3.1 (6.10), (12.5), (12.6) and (12.10) the Lemma follows. 



(12.10) 
□ 



Lemma 12.2 For any B > 0, the function h defined by means of: 



2iZ 
A - 1 



1 + In t - In 



2*Z 



A - 1 



+ ln(2\Are ig t 



lnf|^-iln 



2iZ 
A- 1 



(12.11) 



satisfies 



\h(£,Z,t)\ <C [Z 2 + 



for \Z\ < -^p-, Z £ D(£,B) and |£| > £0 f or some positive constants C = C{B) and £0 = Co(-B) 
sufficiently large. 



Proof. This Lemma is a direct consequence of Lemma (12.1). The only difficulty in order to 
estimate the function $ defined in (6.10) comes from the term V(£)/V(£ + Y) which corresponds to 

y also be written as: 

ln(-$(r}))@(r}-ti,y/\l\Z)dv ■ 



the integral term. This term, given by formula (5.15) may also be written as: 

v(0 



A(£, Z) = — = exp 

v(t+^W\z) 



2 

(A-l)i 



/rrt rj=f3i 



□ 



As for the critical point Z c of <&(£, -,i), the precise result is the following. 

Lemma 12.3 Suppose that B > 2^/tt. Let us define 5q = (A — For all < t < 1 arzd 

Im£ € ((3q, Pq + (A — l)/2), |£| t 2 > £q(B) with £q(B) a positive constant sufficiently large, there 
exists a unique point Z c £ D(£, B) \ Bs t(0) such that <9<J?(£, Z c , t)/dZ = 0. Moreover the following 
asymptotics holds uniformly for < t < 1: 



2iZ c 
A- 1 



2nt{l + iQ)) 1 + 



1 



as \Re£\ t 2 — ► co. 



(12.12) 
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Proof. Computing the derivative of the function $ given by (6.10) gives for all Z £ D(£, B): 
d<$> 



dZ 



(A -1)2 



In (—$(77)) 



C 2 



|f| I + I In ^ i(Z/t) 



1 _ e #rM«lz-i+0 



1 + Am if i 1/2 



jdr] 



X-lJ \ X-l J 2yf\l\Z A-l 



We compute the leading term of the integral in the right hand side of (12.13) as |£| — > 00: 



T (VK|z-»j+£) 



I(£,Z) = / In (-$(??)) 

•^2 A _ e ^(^f^+fA 



dry 



C-2 



In f-*(a + £ + y/\£\Z) 



4tt \ 2 

l-e _ ^ CT 



;i2.i3) 



where C2 = C2 — £—\/|£[.Z. Using Proposition 3.1 we have that, uniformly for Z £ D(£, B), \Z\ < B: 

In + e + vl^)) = In (-$(£ + vlOT) + ^(0 ||| + O (||^) 

for |£| — ► +00, where A(£) is a bounded function of sign(Re^) . It then follows: 



7(£,Z)=ln(-<D(£+v^) 



A(0 



_ 47T 



Ca ( 1 - e ^-1 _X " 



a ^ dfT + L ° ( j^2 



C 2 h_ e -^T' 



c 2 



da + 



_ 4-7T 



4ir \ 2 
1 _ g A^T 



+0 (e~ a|51 J as |e| -> +00 
uniformly for Z £ D(£,B), \Z\ < B, where o > is independent on £, Z. Using 



47T 



_ 47T 

e 



da 



l-e~—i a 



we obtain: 



_ 47T 

e A -i 



c 2 h 



A-l 

47T 



(12.14) 



On the other hand, in order to estimate the second term in the right hand side of (12.14) we deform 



the integration contour C2 to a horizontal line at a bounded distance of the real axis. The resulting 
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integral con then be bou nded by a positive constant independent of £, Z. The third term in the 
right hand side of (12.14) can be bounded using the specific form of C2 as: 



C 



\Z\ 2 \Z\ 3 
+ 



lei lei 1 / 2 



We notice that by Proposition 3.1 we have: 



In (-*(£ + s/W\Z)) = In (-H0) + ol^ 



\z\ 



as |£| — ► +00 



uniformly for Z G D(£,B), \Z\ < B. Combining everything we deduce: 



A-l 
4tt 



In (-$(£)) + 



\Z\ 



ieiva 



as |£| — ► +00 



(12.15) 



uniformly for Z G D(£,B), \Z\ < B. Combining fiTl^ and fllTTo) it follows: 



9$ 



(A-l) 



In (-$(£))+ In 



a-W 4^z v ?l ; 

as iei 



+ 



+0 



|Z| 

W 2 



+00 



uniformly for Z G D(£,B), \Z\ < B. 

Using Rouche's Theorem it then follows that for |£|i 2 sufficiently large, £ G Franja, < t < 1, 
there exists a unique root of (d$>/dZ)(£, Z, i) = in Z G D(£,B) \ B$ ot satisfying the asymtotics 
(12.3) and the Lemma follows. □ 

We now derive estimates for higher order derivatives of <£. 



Lemma 12.4 Suppose that Z c is as in Lemma 12.3. Then the following asymptotics holds: 

1 ^ ' 



(A - 1) Z c 



1 + 



Kit, 



dz 2 

uniformly in < t < 1. 

Proof. The second derivative of $ with respect to Z is: 



as \£\t z 



+00. 



d 2 $ __ 8vr|£|i 
9Z 2 ~ "(A-l) 2 

8vr|£|i 



pn(-*fa))]' 



T (VI«l-'?+0 



2\/W 



c 2 



(A-l) 2 



$'fa) 
c 2 *(»?) 



1 - e< 



(A - 1)2 2Z 2 ' 



Taking Z = Z c and using Lemma 12.3 we deduce, uniformly in < i < 1 



d 2 $ 
SZ 2 



(C,^,i) 



2»viei 



8vr|£|i 



(A-1)Z C (A -1)5 



T^Vl^l-'J+O 



1 _ e ^(ZcVW\-v+0) 



[12.16) 



49 



as |£|i 2 — > +00. Using Proposition 3.1 



AO 



+ / o 



$'(r?) 



i(ZW\t\-V+t) 



-,dr\ 



1 I' 1 



jdrj 



1 



j(Z c y/\i\-vH) 



c 2 \l£l 2 / A _ ^(ZcViCHi-Krt 



cfy + C fe- Q ^l) , as + 



-OO 



(12.17) 



where a is a positive constant independent on £. The first term in the right hand side is estimated 
as: 



1 e x^x(ZcVW\-n+0 



jdr] 



1 



+ / O 



c 2 



l£l 2 



T(Zcy/\t\-nrt) 



(1 - e &v°VM-v+Qy 



dr] + o(e- a ^) as |f| -» +00. 



The first term in the right hand side can be computed explicitly and gives —(A — l)/4ir£. The 
second one is estimated using the form of the contour C2 and is bounded by O (\Z c \/\^\). 



Therefore, the first term in the right hand side of (12.17) is estimated as 0(1/|£|) as |£| — > +00. The 
second term is bounded using similar arguments by 0(1/|£| 3 / 2 ) as |£| — > +00. It then follows that 

J(0 = O (J^j as Id - +00. 



Using (12.16) we obtain 
d 2 $ 



+00 



whence Lemma 112.41 follows. 



□ 



Lemma 12.5 Suppose that Z c , B and 5q are as in Lemma 12.3. Then the following asymptotics 
holds: 



d 3 $ 



dz 3 



o 



t 2 



as |£| t 2 — > +00. 



uniformly in < t < I, Z € D(£, B), \Z\ <B. 
Proof of Lemma 112.51 



d 3 $ 
az 3 



(A - I) 2 Jc 2 V Hv) J ( 1 _ ^(z^i-n+o 



(A - 1) Z 2 Z 3 



[12.18) 
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The last two terms of (12.18) are bounded as C{y/\£\/t 2 + 1/t 3 ). and this can be estimated as 
C\f\£\/t 2 for |£|i 2 >> 1. On the other hand we may bound the first term in the right hand side of 



(12.18) using 



$(t?) 



< 



C 



i + M 



the form of the contour The term under consideration is then bounded by a constant. Combining 
all these estimates for the terms in the right hand side of (12.18) the Lemma follows. □ 
Combining Lemma 12.1 and Lemma 12.3 it follows 

Corollary 12.6 For all < t < 1 : 



1 



\t\V2nt(l + iQ)-^hx[\t\V 2 



In 



2iZ c 
A-l 



+ 0(1) as \(,\t 2 -> +00. 



Proof. Using (12.12) and (12.11) we deduce 



2iZ c 
A-l 



1 + ln(t) - In 



2iZ c 
A-l 



+ ln(2Vvre iQ ^ 



In 



2iZ c 
A-l 



1 



|e|V27rt(l + ^)-^l n (|^|V2 



In 



2iZ c 
A-l 



+ 0(1) 



as ieit 2 



+00. 



+ 0(1) 



;i2.i9) 



□ 



Remark 12.7 We must emphasise that the sign of the real part of the main term in the asymptotic 
expansion of & as t 2 |£| — > 00 is fundamental for the construction of our solutions. 

In the next Lemma we prove that the function Z, t) is "well behaved" in a region \Z — Z c \ > St 
and \Z\ < Mt for M > large. 

Lemma 12.8 For a// 5 > and M > large, there exists Eq > and L > suc/i i/iai the function 
$(£, Z, i) satifies: 



when Z lies in the curve 



Re $(£, Z, t) < i?e *(£, Z c , t) - e t 



7 (M)= 7l (M)U72(M)U73(M)\{Z; \Z - Z c \ <5t} 



where 



7l (M) = {Z 
72 (M) = {Z 
73 (M) = {Z 



Z = Z C + A, AGM, |A| <Mi} 

Z = Z c + Mi + Ai, Ae [0, \ImZ c \ + 71]} 

Z = Z c -Mt + Ai, Ae [0, |JmZ c | + 7 i]} 



/or a// £ and i s«c/i i/iai |£| i 2 > L. 
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Proof. Let us consider the auxiliary function: 



e(£,n,t) 



2V7rv l£l n t [i - in(n) + in(n t 



1 



A - 1 ^/vri 



(12.20) 
(12.21) 

(12.22) 



By Lemma 12.1 



m z, t) = e(e, n, t) - ^ in ( lei 1/2 ) + hz, z, t) . 



Moreover, it is easily checked that Oo is the critical point of the function 0(£, O, t). By Lemma 12.3 
we already know that Z C! the critical point of the function $, converges to Qo as |£| t 2 — > +00. We 
are now going to study the behaviour of the function along the curve obtained from j(M) using 



the change of variable (12.21). Due to the convergence properties of the function $ and its critical 
point Z c when |£| t 2 — > +oo this will be enough in order to prove the statement in Lemma 12.8 
We first consider the curve corresponding to 71 (M). It is then enough to consider Re 0(£, £1, t) along 
the points: = S7q + i, A £ M. A straightforward calculation yields: 



27iVie|t#r) 



a = 1 + Q\ 

ip(a) = 1 - - In (l + a 2 ) + - In 2 - a (J- - arctga^j . 

Since ip'(cr) = — f + arctg(a) the point <7o = 1 is a strict minimum for the function ip. It is also 
easily checked that 



i>{a) 



7T 



as a — » +00 



3vr 



as a 



-00. 



It follows that, for any 5 > there exists e > such that if \a — 1| > 5, 

i>(a)-^{l)>e. 

Arguing by continuity this yields the statement of the Lemma when Z lies in the curve 71 (M). 
The evolution of the function 9 along the two curves corresponding to 72 (M) and 73 (M) is studied 
with very similar arguments. We then only consider the case of the curve 72(M). It is then enough 
to consider i?e0(£, along the points: 



Im (re itp - fi ) 



M 



■tt(A-I) 



a straightforward calculation gives: 

itee(£,fi,t) = -2vWlfl n * 



r cos (p(l — In r) — r sin (/? 
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For Z G 72(M) and the constant M sufficiently large, we have that <p > tt/4 + 5, rcos<p < 2 and 
r > M/2. Similarly, if Z G 73 (M) and the constant M sufficiently large, we have that <p < —tt/4 — 5, 
rcosip < 2 and r > M/2. Therefore, we have in both cases: 



r cos cp(l — In r) — r sin (p 



ttQ 



<p) > e > 0. 



for M large enough. 

Using again the convergence properties of the function and its critical point Z c when |£| t 2 — ► 
+00 this yields the statement in Lemma 12.8 when Z lies in the curves j2(M), js(M). □ 

In the following Lemma we extend the behaviour of Z,t) to the region \Z — Z c \ > 5t and 
\Z\ < £iy/\£\ for some £\ > sufficiently small. 



Lemma 12.9 For all 5 > and M > large, there exists a > 0, S\ > and L > suc/i that the 
function Z, i) satifies: 

Re^,Z,t) < -a^\\Z\, 
for all Z G C\ such that Mt < \Z\ < £iy/\£\ and all £ and t such that |£| t 2 > L. 



Proof. We only need to check that the function defined in the proof of the previous Lemma 
behaves linearly when |0| — * +00 and i?e(0) remains constant. This follows from 



22e(e(£,n,*)) < 



7T 



3/2 



\Z\t\n\ 



uniformly for £1 = i \Q\ + 0(1). Using (12.21) we deduce 



Re(@(t,n,t))<-^\(i\ x | 



for all Z G Ci such that |Z| > Mt assuming that M > is sufficiently large. Using Lemma 12.1 
have: 



wc 



Re Z, t) < Re(G(£, Sl,t)) + C Z 2 + O 



<-^M\\z\ 



7T 



2(A-1) 



for all Z G Ci such that Mt < \Z\ < £o^/\^\. The result follows for e\ small enoug h and \i\ 1/2 t -* 
+00. □ 



Lemma 12.10 For all B > i/iere exists £0 C > stic/i t/iai 



9^ 



<cj|j±, £ = 1,2 



(12.23) 



/or y = Z^ ; |y| < f , Z G D(£, B) and |i?e(£)| > 



53 



Proof. Differentiating the function 

2 



1) * Sir, 



ln(-$(r)))Q(r)-t,Y)dri 



(A — i) b j Imri=l3l 

2iY , _ 2iY / 2iY 1 



A-l 



ln(t) 



+ 



A-l V A - 1 



In 



2iY 
A-l 



with respect to £ we obtain: then, 



«9* 



— ^— / In (-$(,7)) ~(v~^ Y)drj 
(A - 1) i Ji m i t= fa d£ 

-JY^V- I h H(l))?(l-^)<i>l 

(A- l)t Jimr,^ Or] 



By Proposition 3.1 we have 



(A — 1) « Jl mr) =f3i HV) 
*'(»?) . c 



*(»?) 



< 



1 + 



in the domain Y/y |£| G D(£,B). We deform the contour of integration to C 2 defined in Figure 4. 
Then we split the integral in two pieces: 



+c 



{i,y,t)\<c 



1 



«7ee a ,|i7-£|<J£l 1 + 



■\Q^-^Y)\\dq\ 



1 



'r?GC2,|r,-CI>^ 1 + 

By the exponential decay of the function 0: 

J 2 < Ce~ a ^ 

for some positive constant a. On the other hand, 

1 



\@{r,-^Y)\\dr,\ = J 1 + J 2 . 



l- e -&*-$ 1 _ e &(-(v-Q+Y) 



The integral J\ is then divided in two parts. The first, J\ t x is the integral along the "vertical part 
of the curve" C 2 . The second, Ji >2 is along the horizontal part of that curve, where Im{rj) = (3\. In 
the integral J^i, Re(rj) is bounded and therefore \@(rj — £, Y)\ < C. 

Since the total length of the integration curve of J%\ is of order \Y\ we deduce that J\ \ < 
C\Y\/(1 + |£|). We split the integral Ji >2 as follows: 



Jl,2 < 



c 



1 + 



Im{r,)=Px,\n-ti\<2\Y\ 



\Q( T] -^Y)\di 1 + I 

Jr. 



Im(v)=Pi,\ri-Z\>2\Y\ 



tY)\drj 



7m(r ? )=ft,|7?-?|>2|r| 



|e(r/-e,y)|^< / 

J I 



/m(j7)=^i,k|>2|r| 



e 



47T 

e 



dry. 
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We use now that, if Re(a) > 2 \Y\: 



47T 



(i 



2 7T | | 



and if i?e(a) < -2|Y| 



e A -! 



47T 

(1 - e 



< C 



r )(l-e- 

The last remaining term is easily estimated by: 



4tt y- 



47T _ 47T /_ V"\ 

(e^ ff - l)(e^ (<T ~ y) - 1) 

4tt y Gtt I _ I 



/ |0(77-^,^)||rf77| < C / 



|dr7| < c|y|. 



/Jm(f7)=/3i,[»7-$|<2|y| 

It then follows that J\ 2 < Ce _a '^'/(1 + |£|) for positive constant C and a. This ends the proof of 



fll2.23[ ) for ^ = 1. Similarly. 

d 2 ^ 



d£ 2 



with 



(A-l) 



e( v -t,Y)d v 



for y/\/|£| € D(£,B), again by Proposition 



< 



C 



arguments as those of the proof of (12.23) for t = 2. 



3.1 



The proof of (12.23) follows then from the same 

□ 
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